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Abstract 

This paper proposes a formal model for a network of robotic agents that move and communicate. 
Building on concepts from distributed computation, robotics and control theory, we define notions 
of robotic network, control and communication law, coordination task, and time and communication 
complexity. We then analyze a number of basic motion coordination algorithms such as pursuit, 
rendezvous and deployment. 

1 Introduction 

Problem motivation The study of networked mobile systems presents new challenges that lie at 
the confluence of communication, computing, and control. In this paper we consider the problem of 
designing joint communication protocols and control algorithms for groups of agents with controlled 
mobility. For such groups of agents we define the notion of communication and control law by extending 
the classic notion of distributed algorithm in synchronous networks. Decentralized control strategies 
are appealing for networks of robots because they can be scalable and they provide robustness to 
vehicle and communication failures. 

One of our key objectives is to develop a computable theory of time and communication complexity 
for motion coordination algorithms. Hopefully, our formal model will be suitable to analyze objectively 
the performance of various coordination algorithms. It is our contention that such a theory is required 
to assess the complex trade-offs between computation, communication and motion control or, in other 
words, to establish what algorithms are scalable and practically implementable in large networks of 
mobile autonomous agents. The need for modern models of computation in wireless and sensor network 
applications is discussed in the well-known report 

Literature review To study complexity of motion coordination, our starting points are the standard 
notions of synchronous and asynchronous networks in distributed and parallel computation, e.g., see 
Lynch and, with an emphasis on numerical methods, Bertsekas and Tsitsiklis j^. This established 
body of knowledge, however, is not applicable to the robotic network setting because of the agents' 
mobility and the ensuing dynamic communication topology. 
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An important contribution towards a network model of mobile interacting robots is introduced 
by Suzuki and Yamashita |4^, see also P, [i^]. The Suzuki- Yamashita model consists of a group 
of "distributed anonymous mobile robots" that interact by sensing each other's relative position. 
A related model is presented in [l^ . A brief survey of models, algorithms, and the need for 

appropriate complexity notions is presented in [4o| . This model is discussed in Section [2.41 

Recently, a notion of communication complexity for control and communication algorithms in 
multi-robot systems is analyzed in see also (2ol | where a formal model of communication and 



control laws for multi-agent networks is proposed. A general modeling paradigm is discussed in 2^|. 
The time complexity of a distributed algorithm for coordinated motion planning is computed in 41 1 . 

Finally, we conclude this literature review session with a general overview of the recent large 
research effort dedicated to multi-vehicle systems in the robotics and control domains. A survey 
on cooperative mobile robotics is presented in 0] and an overview of control and communication 
issues is contained in 21 1 . Specific topics of interest include formation control 0, 0, ll^. jj^ l43 ] , 
rendezvous 0, 0, 2L , , flocking , swarm aggregation 0] , gradient climbing |3,'j | , cyclic 
pursuit li^, robotic exploration [i^, deployment 0,a, and forag^g HEiConsensus and control 



iheoretic^problems on dynamic ^aphs are discussed in UW^-d^, respectively. 



Statement of contributions We summarize our approach as follows. A robotic network is a group 
of robotic agents moving in space and endowed with communication capabilities. The agents position 
obey a differential equation and the communication topology is a function of the agents' relative po- 
sitions. Each agent repeatedly performs communication, computation and physical motion as follows. 
At predetermined time instants, the agents exchange information along the communication graph and 
update their internal state. Between successive communication instants, the agents move according to 
a motion control law, computed as a function of the agent location and of the available information 
gathered through communication with other agents. In short, a control and communication law for 
a robotic network consists of a message-generation function (what do the agents communicate?), a 
state-transition function (how do the agents update their internal state with the received informa- 
tion?), and a motion control law (how do the agents move between communication rounds?). We 
then define the notion of time complexity of a control and communication law (aimed at solving a 
given coordination task) as the minimum number of communication rounds required by the agents to 
achieve the task. The time complexity of a coordination task, as opposed to that of an algorithm, is the 
minimum time complexity of any algorithm achieving the task. We also provide similar definitions for 
notions of mean and total communication complexity. In particular, we show that our model of time 
and communication complexity satisfies a basic well-posedness property that we refer to as "invariance 
under reschedulings." 

Building on these notions we establish complexity estimates for a few basic motion coordination 
algorithms such as pursuit, rendezvous, and deployment. First, for a network of agents moving on 
the circle, we introduce a dynamic law, called agree-and-pursue, that combines elements of the classic 
cyclic pursuit and leader election problems; see |27l. Eil. respectively. We show that this law achieves 
consensus on the agents' direction of motion and equidistance between the agents' positions. Further- 
more, we show that these tasks are complete in time 0{N) and 0{N'^ log A^), respectively. Second, we 
analyze a simple averaging law for a network of locally-connected agents moving on a line. This law is 
related to the widely known Vicsek's model, see We show that this law achieves rendezvous 

(without preserving connectivity) and that its time complexity belongs to Vt{N) and 0{N^). Third, 
for a network of locally-connected agents moving on a line or on a segment, we show that the well- 
known circumcenter algorithm by []| has time complexity of order @{N). (This algorithm achieves 
rendezvous while preserving connectivity with a communication graph with 0{N'^) links.) We then 
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consider a network based on a different communication graph, called the limited Delaunay graph, that 
arises naturally in computational geometry and in the study of wireless communication topologies. 
For this less dense graph with 0{N) communication links, we show that the time complexity of the 
circumcenter algorithm grows to @{N'^ \ogN). For a network of agents moving on M"^ (with a certain 
communication graph) we introduce a novel "parallel-circumcenter algorithm" and establish its time 
complexity of order Q{N). Fourth and last, for a network of agents in a one-dimensional environment, 
we show that the time complexity of the deployment algorithm introduced in is 0{N^ log N). To 
obtain these complexity estimates, we develop some novel analysis methods. In particular, we develop 
a key set of results on tridiagonal Toeplitz and circulant matrices that characterize their convergence 
rates as a function of the matrices dimensions. 

We refer the reader to the Conclusions (Section [SJ for a discussion about the limitations of the 
proposed results and about avenues for future research. 

Organization Section |2l presents a general approach to the modeling of robotic networks by for- 
mally introducing various notions including, for example, those of communication graph, control and 
communication law, and network evolution. Section 121 defines the notions of task and of time and 
communication complexity for a control and communication law. Section 0] presents a few motion co- 
ordination algorithms performing the basic tasks of rendezvous and deployment. For each algorithm, 
we characterize time and communication complexity, along with asymptotic correctness. Finally, we 
present our conclusions in Section [51 In the appendices, we review some basic computational geo- 
metric structures and we prove some key facts about the discrete-time dynamical systems defined by 
tridiagonal Toeplitz and circulant matrices. 

Notation We let BooleSet be the set {true, false}. We let rTie{i n] '^i denote the Cartesian 
product of sets Si, ... , Sn- We let M+ and M+ denote the set of strictly positive and non-negative real 
numbers, respectively. We let N denote the set of positive natural numbers and No denote the set of 
non-negative integers. If iS is a set, then diag(5 x S) = {{s,s) £ S x S \ s G S}. For x € M, we let 
[x\ denote the floor of x. For x G M"', we let ||x||2 and || x||qo denote the Euclidean and the cx)-norm 
of X, respectively. Recall that ||x||oo < ||a;||2 < V^||2;||oo for all x G W^. For x € M*^ and r G M+, we 
let B{x,r) and B{x,r) denote the open and closed ball in centered at x of radius r, respectively. 
We let ei, . . . , be the standard orthonormal basis of W^. Also, we define the vectors = (0, . . . , 0) 
and 1 = (1,... ,1) in R'^. For f,g: N ^ M, we say that / G 0{g) (respectively, / G ^{g)) if there exist 
A^o G N and A: G R+ such that |/(iV)| < k\g{N)\ for ah N > Nq (respectively, |/(A^)| > k\g{N)\ for all 
A^ ^ A^q). If / G 0{g) and / G ^{g), then we use the notation / G Q{g). Finally, we refer the reader 
to Appendix El for some useful geometric concepts. 

2 A formal model for synchronous robotic networks 

In this section we introduce a notion of robotic network, as a group of robotic agents with the ability 
to move and to communicate according to a specified communication topology. Our notion of control 
and communication law for a robotic network parallels the definition in |27| of synchronous distributed 
algorithm. 

2.1 The physical components of a robotic network 

Here we introduce our basic definition of physical quantities such as the agents and such as the ability 
of agents to communicate. We begin by defining a control system as our basic model for how each 
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robotic agent moves in space. 

Definition 2.1 A control system is a tuple {X,U, Xq, f) consisting of 
(i) X is a differentiable manifold, called the state space; 
(a) U is a compact subset o/M™ containing 0, called the input space; 
(Hi) Xq is a subset of X, called the set of allowable initial states; 

(iv) f: X X U ^ TX is a C^-map with f{x,u) G T^X for all (x,n) £ X xU. • 

We refer to x £ X and u £ U as a state and an input of the control system, respectively. We will 
often consider control-affine systems, i.e., control systems for which f{x,u) = /o(x) + X^^i fa{x)ua- 
In such a case, with a slight abuse of notation, we will represent the map / as the ordered family of 
C°°-vector fields (/o, /i, . . . , fm) on X. We will also sometime consider driftless systems, i.e., control 
systems for which /(x,0) = 0. 

Definition 2.2 A network of robotic agents (or robotic network^ S is a tuple (/, E'cmm) consisting 
of 

(i) / = {!,... , A^}; / is called the set of unique identifiers (UIDs); 

(a) A = {^Wjig/ = {(XW, [/['I, Xot*^, /f*^)}iG/ is a set of control systems; this set is called the set of 
physical agents; 

(Hi) Ecrara is a map from Jlie/"'^^*^ ^'^ subsets of I x I \ diag(/ x /); this map is called the 
communication edge map. 

If A^^^ = {X,U, Xq, f) for all i £ I, then the robotic network is called uniform. • 

Let us comment on this definition and on how robotic agents communicate in a robotic network 

{,1 J A, -E'cmm) • 

Remarks 2.3 (i) By convention, we let the superscript [i] denote the variables and spaces which 
correspond to the agent with unique identifier i; for instance, xW G and Xq^ € X^^ denote 
the state and the initial state of agent ^W, respectively. We refer to {x^^\ . . . ,x^^^) G riig/^^*^ 
as a state of the network. 

(ii) The map -Ecmm models the topology of the communication service between the agents. In 
other words, at a network state x = {x^^\ . . . , x^^^) , two agents at locations and x^'^ can 
communicate if the pair {i,j) is an edge in Ecmm{x^^h ■ ■ ■ , x^^^). Accordingly, we refer to the pair 
(/, Ecmm{x^^h • • • ) x^^^)) as the communication graph at x. When and what agents communicate 
is discussed below. 

Maps of the form E: Yli^j X^^^ 2^^^\^^^sii^i) ^re called proximity edge maps, they arise in 
wireless communication and in computational geometry, we refer the reader to Appendix 1X1 for 
more details. Excluding edges of the form {i,i), for i G I, means that any individual agent does 
not communicate with itself. • 

To make things concrete, let us now present some useful examples of robotic networks; we will 
use these examples throughout the paper. We start with a fairly common example and define some 
interesting variations. 
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Example 2.4 (Locally-connected first-order agents in M^) Consider N points x'^^ , . . . , x^^l in 

the Euclidean space M"', d > 1, obeying a first-order dynamics = uH(t). According to Defi- 

nition [23 these are identical agents of the form A = (M'^, M"', M'^, (0, ei, . . . ,ed)). Assume that each 
agent can communicate to any other agent within Euclidean distance r, that is, adopt as communi- 
cation edge map the r-disk proximity edge map -Er-disk defined in Appendix These data define the 
uniform robotic network S-^d j._(jisk ~ 

-disk ) ■ • 

Example 2.5 (LD-connected first-order agents in M'^) Consider the set of physical agents de- 
fined in the previous example. For r G R+, recall from Appendix 1^ the r-limited Delaunay map -Br-LD 
defined by 

£ Er-LDix^^K- ■ • if and only if {V^^nB{x^^, §)) n {V^^^ nB{x^^\ ^)) / 0, i / j, 

where {T^^^', . . . , yt^l} is the Voronoi partition of M'^ generated by {x'^^\ . . . , x'^^}. These data define 
the uniform robotic network S^d = (Ij^A., -Br-Lo)- • 

Example 2.6 (Locally-oo-connected first-order agents in M'^) Consider the set of physical agents 
defined in the previous two examples. For r G R+, define the proximity edge map -Er-oo-disk by 

G -B,._oo-disk(2;W, . . • if and only if ||x[*l - x[^'l||oo < r, i / j. 

These data define the uniform robotic network S^d ^.oo-disk ~ (^)-^i -E'r_oo-disk)- • 

In what follows, let be the unit circle, and measure positions on the circle counterclockwise from 
the positive horizontal axis. Let us be specific about distances on the circle and related concepts. For 
x,y G we let dist(x,y) = min{distc(x, y), distcc(ic, y)}- Here, distc(x,y) = {x — y) (mod27r) is the 
clockwise distance, that is, the path length from x to y traveling clockwise. Similarly, distcc(2;,y) = 
{y — x) (mod27r) is the counterclockwise distance, i.e., the path length from x to y traveling counter- 
clockwise. Here x (mod27r) is the remainder of the division of x by 2tt. 

Example 2.7 (Locally-connected first-order agents on the circle) For r G M+, consider the 
uniform robotic network 5gi ^-disk = (-^i Er-disk) composed of identical agents of the form (S^, (0, e)). 
Here e is the vector field on describing unit-speed counterclockwise rotation. We define the r-disk 
proximity edge map -Br-disk on the circle by 

G ^r-disk(^^^^ ■ ■ • if and only if dist{9^\ 9^^^) < r , 

where dist(x, y) is the geodesic distance between the two points x, y on the circle. • 

2.2 Control and communication laws for robotic networks 

Here we present a discrete-time communication, continuous-time motion model for the evolution of a 
robotic network. In our model, the robotic agents evolve in the physical domain in continuous-time 
and have the ability to exchange information (position and/or dynamic variables) that affect their 
motion at discrete-time instants. 

Definition 2.8 (Control and communication law) Let S be a robotic network. A (synchronous, 
dynamic, feedback) control and communication law CC for S consists of the sets: 

(i) T = {t^l^GNo C M4. is an increasing sequence of time instants, called communication schedule; 
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(ii) L is a set containing the null element, called the communication language; elements of L are 
called messages; 

(Hi) Vl^W, i £ I, are sets of values of some logic variables w'*', i G I; 

(iv) VFg*^ C W^^\ i I, are subsets of allowable initial values; 
and of the maps: 

(i) msgf*^ : T x x x I ^ L, i £ I , are called message-generation functions; 

(ii) stfl^l : T X W^W X ^ W^^, i G /, are called state-transition functions; 

(Hi) ctlW : S+ X X X W^^ x ^ U^^, i G /, are called control functions. • 

We will sometimes refer to a control and communication law as a motion coordination algorithm. 
Control and communication laws might have various properties. 

Definition 2.9 (Properties of control and communication laws) Let S be a robotic network 
and CC be a control and communication law for S. 

(i) If S is uniform and if W^'^^ = W, msgW = msg, stfW = stf, ctl^^^ = ctl, for all i £ I, then CC is 
said to be uniform and is described by a tuple (T, L, W^, msg, stf, ctl). 

(ii) If = Wq^ = for all i £ I, then CC is said to be static and is described by a tuple 
(T,L,{msgH}ie/,{ctl[*l}ie/), mi/i msgW : TxX^^xI ^ L, and ctl'*! : TxX^ xXW xL^ ^ f/W. 

(Hi) CC is said to be time-independent if the message-generation, state-transition and control func- 
tions are of the form msgW : x W^^ x I ^ L, stf'^l : W^^ x ^ W^^ , ctl^ : x X^ x 
Vl^H X C/H, i G /, respectively. • 

Roughly speaking this definition has the following meaning: for all i £ I, to the ith physical agent 
corresponds a logic process, labeled i, that performs the following actions. First, at each time instant 
ti G T, the ith logic process sends to each of its neighbors in the communication graph a message 
(possibly the null message) computed by applying the message-generation function to the current 
values of and w^^h After a negligible period of time (therefore, still at time instant t^ G T), the 
ith logic process resets the value of its logic variables by applying the state-transition function to 
the current value of w'*', and to the messages received at time t£. Between communication instants, 
i.e., for t G [t£, t^+i), the motion of the ith agent is determined by applying the control function to the 
current value of x'*', the value of a;W at t£, and the current value of w^^\ This idea is formalized as 
follows. 

Definition 2.10 (Evolution of a robotic network subject to a control and communication law) 

Let S be a robotic network and CC be a control and communication law for S. The evolution of {S,CC) 
from initial conditions Xg' G Xq'*' and Wq^ G Wq\ i £ I, is the set of curves x^'^: XW, 
i G /, ^ G No, and w^^ : T M^H, i £ I, satisfying 

i^U{t) = /(xW'^(t), ctlW(t,xH'^(t),xH'^(t,),«;H(t,),yH(t,))), 

where, for £ G No, and i £ I, 

xW'^(t,) = xH'^-i(t,) , u;H(t,) = stf[^l(t,,i^;H(t,_i),y[*l(t,)) , 



6 



with the conventions that x^' ^(to) = a^o' and = 'Wq\ i ^ I. Here, the function ; Y 

(describing the messages received by agent i) has components yf{ti), for j G /, given by 



L 



N 



\msg^^\t,,xl^^''-\t,),wl^\t,_,),i), zf {i,j) G E,^^{xm''-Ht,),...,xW''-\t,)), 
I null, otherwise. 



Let us emphasize two limitations regarding the proposed communication model. 

Remarks 2.11 (Idealized aspects of communication model) (i) We refer to CC as a syn- 
chronous control and communication law, because the communications between all agents takes 
place always at the same time for all agents. We do not discuss here the important setting of 
asynchronous laws. 

(ii) The set L is used to exchange information between two robotic agents. The message null 
indicates no communication. We assume that the messages in the communication language 
L allow us to encode logical expressions such as true and false, integers, and real numbers. 
A realistic assumption on the communication language would be to adopt a finite-precision 
representation for integers and real numbers in the messages. Instead, in what follows, we 
neglect any inaccuracies due to quantization and we discuss this topic as an open problem in 
Section |3 • 

In the following remarks we introduce additional useful notations and emphasize an important 
assumption. 

Remarks 2.12 (Related concepts and notations) (i) To distinguish between the null and the 
non-null messages received by an agent at a given time instant, it is convenient to define the 
natural projection ttl : 2^ that maps an array of messages y to the subset of L containing 

only the non-null messages in y. 

(ii) In many uniform control and communication laws, the messages interchanged among the network 
agents are (quantized representations of) the agents' states and dynamic states. The correspond- 
ing communication language is L = X x W and message generation function msgg^^j : T x X x 
WxI^XxW is referred to as the standard message-generation function and is defined by 
msgstd(*>2;,'u;, j) = {x,w). 

(iii) By concatenating the curves x^'^ and tti^'^, for ^ E No, we can define the evolution of the ith 
robotic agent M.+ 3 t ^ {x^'\t) , w^^ (t)) e X^^ x W^^. Additionally we can define the curves 

W+Bt^x{t) = (xt^l (t) , . . . , xl^l (t)) G n ' 

is/ 

W+3t^w{t) = (w;W(t),...,«;[^l(t)) G JJl^W. 

(iv) In the proposed notion of synchronous robotic network, we assume that the states (xl^l , . . . , x^^l) G 
J^jgj-XW are not reset after a communication round, i.e., we are setting xW'^(t^) = x^^^'^~^{ti). 
This is reasonable for instance when (x'^^ , . . . , x^^l) are positions or other physical variables. It 
would be possible, though, to allow for more general situations by defining state-transition func- 
tions of the form stf^^l : T x X^^ x W^^ xL^ ^ x W^^, for i G /, and setting (x^-^ w^^^^){ti) = 
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st&\ti,x^^^'^~^{ti),w^^^'^~^{ti),y^'^^{ti)). Other possible extensions would be to consider (1) non- 
deterministic continuous-time evolutions of the state variables (by means of controlled differential 
inclusions, as opposed to the controlled differential equations) or (2) non-deterministic discrete- 
time updates of the dynamical variables (by means of set- valued state-transitions functions). 
Further generalizations are also possible, but we do not consider them here in the sake of sim- 
plicity. • 



2.3 Example control and communication laws 

It is now a good time to present various examples of control and communication laws. We start by 
considering a very simple averaging algorithm as a static control and communication law. We then 
present an interesting dynamic control and communication law on the circle. 

Example 2.13 (The move-toward-average control and communication law) From ExampleEU 
consider the uniform network S^d j._disk °^ locally-connected first-order agents in W^. We now define a 
static, uniform and time-independent law that we refer to as the move-toward-average law and that 
we denote by CCavrg- We loosely describe it as follows: 

[Informal description] Communication rounds take place at each natural instant of time. 
At each communication round each agent transmits its position. Between communication 
rounds, each agent moves towards and reaches the point that is the average of its neighbors' 
positions; the average point is computed including the agent's own position. 

Next, we formally define the law as follows. First, we take T = No and we assume that each agent oper- 
ates with the standard message-generation function, i.e., we set L = W^ and msg{x,j) = ms,g^^^{x , j) = 
X. Second, we define the control function ctl: 

ctl(x, Xsmpid, y) = -kprop vers {x - avrg(y U {xsmpid})) , 

where fcprop > r, vers: M'^ — > M'^ is defined by vers(O) = and vers(f) = w/||f||2 for v ^ 0, and the 
map avrg computes the average of a finite point set in M'^: 

avrg(5) = — ^ p . 

1 p£n^(S) 

In summary we set CCavrg = (No, I^"^, nisgs^di ctl). An implementation of this control and communica- 
tion law is shown in Fig. ^ for d = 1. Note that, along the evolution, (1) several agents rendezvous, 
i.e., agree upon a common location, and (2) some agents are connected at the simulation's beginning 
and not connected at the simulation's end. Finally, we remark that this law is related to the Vicsek's 



model discussed in |l7l . |44 1 . 



Example 2.14 (The agree-and-pursue dynamic control and communication law) From Ex- 
ample 12.71 consider the uniform network cSgi ^.jisk of locally-connected first-order agents in . We 
now define the agree-and-pursue law, denoted by CCagr-pursuit, as the uniform and time-independent 
law loosely described as follows: 

[Informal description] The dynamic variables are drctn taking values in {c, cc} and prior 
taking values in /. At each communication round, each agent transmits its position and 
its dynamic variables and sets its dynamic variables to those of the incoming message 
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Figure 1: Evolution of a robotic network under the move-toward-average control and communication 
law in Example 12. LSI implemented over the r-disk graph, with r = 1.5. The vertical axis corresponds 
to the elapsed time, and the horizontal axis to the positions of the agents in the real line. The 51 
agents are initially randomly deployed over the interval [—15, 15]. 



with the largest value of prior. Between communication rounds, each agent moves in the 
counterclockwise or clockwise direction depending on whether its dynamic variable drctn 
is cc or c. For /cprop G]0, ^[, each agent moves /cprop times the distance to the immediately 
next neighbor in the chosen direction, or, if no neighbors are detected, fcprop times the 
communication range r. 

Next, we formally define the law as follows. Each agent has logic variables w = (drctn, prior), where 
wi = drctn G {cc, c}, with arbitrary initial value, and W2 = prior E I, with initial value equal to 
the agent's identifier i. In other words, we define W = {cc, c} x /, and we set Wq^ = {cc, c} x {i}. 
Each agent i G I operates with the standard message-generation function, i.e., we set L = §1 X and 
msgW = msggtj, where m.sgs^^^{9 , w , j) = {9,w). The state-transition function is defined by 

sti{w,y) = argmax{z2 | z G {iTL{y))2 U {w}]. 
For /cprop G the control function is 



ctl(6',6'smpid,u^,y) = Vop 



min{r} U {distcc(6'smpld, 6'rcvd) I ^'rcvd G (vrL(y))i}, if drctn = cc, 
-min{r} U {distc(6'smpld, ^'rcvd) I ^'rcvd G (7rL(y))i}, if drctn = c. 



Finally, we sketch the control and communication in equivalent pseudocode language. This is 
possible for this example, and necessary for more complicated ones. For example, the state-transition 
function is written as: 

function stf ((drctn, prior), y) 

for each non-null message (^rcvdj (drctn^vd) P^'iorj.^^^)) in y: 
if (priorj,pyj > prior), then 
drctn := drctni-cvd 
prior := prior^.^^^ 
endif 
endf or 

return (drctn, prior) 
Similarly, the control function ctl is written as: 
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function ctl (6*, ^smpld, (drctn, prior), y) 

for each non-null message (^rcvd, (drctUfcvd, pi"ioi"j.j,yjj)) in y: 
if (drctn = cc) AND (distcc(6'smpld, ^'rcvd) < dtmp), then 

'^tmp • — distcc (^smpldj ^rcvd) 

else if (drctn = c) AND (distc(6'smpid, ^rcvd) < ^tmp), then 

dtrap '■= distc(^smpld, ^rcvd) 

endif 
endf or 

if (drctn = cc), then return A;propdtmp) else return — A;propdtmp 
endif 

An implementation of this control and communication law is shown in Fig. |21 Note that, along the 
evolution, all agents agree upon a common direction of motion and, after suitable time, they reach 




Figure 2: The agree-and-pursue control and communication law in Example 12. 141 with = 45, r = 
27r/40, and kp^p = 1/4. Disks and circles correspond to agents moving counterclockwise and clockwise, 
respectively. The initial positions and the initial directions of motion are randomly generated. The 
five pictures depict the network state at times 0, 12, 37, 100, 400. 

2.4 Groups of robotic agents with relative-position sensing 

In this last subsection on modeling, we discuss in some detail the Suzuki- Yamashita model mentioned 
in the Introduction, see |4^. This model consists of a group of identical mobile robots characterized 
as follows: no explicit communication takes place between the agents, at each instant of an "activation 
schedule," each robot measure the relative position of all other robots and moves according to a spec- 
ified algorithm. In this model, robots are referred to as "anonymous" and "oblivious" in precisely the 
same way in which we defined control and communication laws to be uniform and static, respectively. 

As compared with our notion of robotic network, the Suzuki- Yamashita model is more general in 
that the robots' activations schedules do not necessarily coincide (i.e., this model is asynchronous), 
and at the same time it is less general in that (1) robots cannot communicate any information other 
than their respective positions, and (2) each robot observes every other robot's position (i.e., the 
complete communication graph is adopted; this limitation is not present for example in 0]). Note 
that a control and communication law, as in our definition, can be implemented on a synchronous 
Suzuki- Yamashita model if the law (1) is static and uniform, (2) only relies on communicating the 
agents' positions (e.g., the message-generation function is the standard one), and (4) entails a control 
function that only depends on relative positions (as opposed to absolute positions). 
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3 Coordination tasks and complexity measures 



In this section we introduce concepts and tools useful to analyze a communication and control law. 
We address the following issues: What is a coordination task for a robotic network? When does a 
control and communication law achieve a task? And with what time and communication complexity? 

3.1 Coordination tasks 

Our first analysis step is to characterize the correctness properties of a communication and control 
law. We do so by defining the notion of task and of task achievement by a robotic network. 

Definition 3.1 (Coordination task) Let S be a robotic network and let W be a set. 
(i) A coordination task for S is a map T: Oie/^'*' X ^ BooleSet. 

(ii) IfW = 0, then the coordination task is said to be static and is described by a map T: Hie/ -^^'^'^ ~^ 
BooleSet. 

Additionally, let CC a control and communication law for S. 

(i) The law CC is compatible with the task T: Hie/"^^'^ ^ ~^ BooleSet if its logic variables 
take values in W, that is, if W^^'^ = W, for all i £ I. 

(ii) The law CC achieves the task T if it is compatible with it and if, for all initial conditions Xq' G Xq^ 
and Wq'^ G ^d''> i £ I, the corresponding network evolution t {x{t),w{t)) has the property that 
there exists T G ]R_|_ such that T{x{t),w{t)) = true for all t>T. • 

Loosely speaking, achieving a task might mean obtaining a specified pattern in the position of the 
agents or of their dynamic variables. We now give some examples of interesting coordination tasks. 

Example 3.2 (Rendezvous tasks) First, let S = {I,A,Ecmm) be a uniform robotic network. The 
(exact) rendezvous task 'T^ndzvs '■ X^ — > BooleSet for S is the static task defined by 



'^ndzvs {x^ \ . . . , X^ ^) 



I true, if = xb'l, for all G Scmm(xW, . . . ,x[^l; 
1 false, otherwise. 



Second, let 5 = (I,^, £^cmm) be a uniform robotic network with agents' state space X C M*^. Ex- 
amples networks of this form are S^d ^.^isk) see Examples 12 . 41 and 12.131 and S^d ^.ldj see Examples E 
For e > 0, the e-rendezvous task T^-mdzvs : X^ — > BooleSet for S is defined by 



'^-rndzvs (^) 



true, if 



- avrg ({xl'l} U {xt-^'l | (i, j) G ^cmm(x)}) 



< e, for all i G /, 

2 



false, otherwise. 



where x = (xW, . . . , xl^l) G C (M'^)^. In other words, T^-mdzvs is true at x G {W^)^ if, for 
all i £ I, xW is at distance less than e from the average of its own position with the position of its 
-E'cmm-neighbors. • 
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Example 3.3 (Agreement and equidistance tasks) From Example 12. 7( consider the uniform net- 
work 5gi j._disk °f locally-connected first-order agents in From Example 12.141 recall the agree-and- 
pursue control and communication law CCagr-pursuit with dynamic variables taking values in W = 
{cc, c} X /. There are two tasks of interest. First, we define the agreement task Tdrctn : (S^)^ x W'^ — > 
BooleSet by 



'7drctn(6','f^) 



{true, if drctnW = . . . = drctnt^l, 
false, otherwise, 



where 6 = (6*^, . . . ,0^^^), w = . . . and = (drctnW,priorW), for i G /. Furthermore, 

for e > 0, we define the static e- equidistance task T^qdstnc : (S^)^ BooleSet by 

_ Jtrue, if I minj^, distc(0[*l, ^I^l) - minj^, distcc{0^^ ,9^^^ < e, for all i G /, 
I false, otherwise. 

In other words, T^-eqdstnc is true when, for every agent, the clockwise distance to the closest clockwise 
neighbor and the counterclockwise distance to the closest counterclockwise neighbor are approximately 
equal. • 

Example 3.4 (Deployment tasks) By optimal deployment on the convex simple polytope Q C 
M.'^ with density function cj): Q ^ ]R_|_, we mean the following objective: place the agents on Q so 
that the expected square Euclidean distance from any point in Q to one of the agents is minimized. 
To define this task formally, let us review some known preliminary notions; we will require some 
computational geometric notions from Appendix EI We consider the following network objective 
function ?^dcplmnt : 

Wdeplmnt(xW,...,x[^])= / min ||g - ||i </)(g)dg . (1) 

This function and variations of it are studied in the facility location and resource allocation research 
literature; see convenient to study a generalization of this function. For r € M+, define 

the saturation function sat^ : M — > M by sa,tr{x) = x if x < r and satr(a;) = r otherwise. For r G M_|-, 
define the new objective function Wr-dcpimnt • 

^ M by 

7^r-depimnt(a;W,...,x[^]) = / min sat . (||g - </,(g)dg . (2) 

Jq 

Note that if r > 2diam((5), then Hdepimnt = "^r-depimnt • Let {F^, . . . ,1/^^!} be the Voronoi partition 
of Q associated with {x'^l, . . . The partial derivative of the cost function takes the following 

meaningful form 

^l^I^^{x^^\. . . = 2Mass(yW nB(xH, §))( Centroid(yW nB(x[^ §)) - x^) , i G J. 

(Here, as in Appendix IXl Mass(S') and Centroid(5) are, respectively, the mass and the centroid of S C 
M"^.) Clearly, the critical points of "H^.dcpimnt are network states where = Centroid(yW n i?(xW, §)). 
We call such configurations |-centroidal Voronoi configurations. For r > 2diam((5), they coincide with 
the standard centroidal Voronoi configurations on Q. Fig. illustrates these notions. 
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Figure 3: Centroidal and |-centroidal Voronoi configurations. The density function (p is depicted by 
a contour plot. For each agent i, the set ni?(pj, |) is plotted in light gray. 



Motivated by these observations, we define the following deployment task. For r, e G M+, define 
the e-r- deployment task T^-r-dcpimnt • 

qN _^ BooleSet by 

Jtrue, if \\x^^ - Centroid(yH n:B(xH, §))||2 < e, for alH G / , 
I false, otherwise. 

Roughly speaking, T^-rdopimnt is true for those network configurations where each agent is sufficiently 
close to the centroid of an appropriate region T^W niJ(xW, |). Note that, with a natural definition of 
Voronoi partitions of the circle, the two tasks 7^_eqdstnc and '7^-r-dcpimnt are closely related. • 

3.2 Complexity notions for control and communication laws and for coordination 
tasks 

We are finally ready to define the key notions of time and communication complexity. These notions 
describe the cost that a certain control and communication law incurs while completing a certain 
coordination task. We also define the complexity of a task to be the infimum of the costs incurred by 
all laws that achieve that task. 

First we define the time complexity of an achievable task as the minimum number of communication 
rounds needed by the agents to achieve the task T. 

Definition 3.5 (Time complexity) Let S be a robotic network and let T be a coordination task for 
S. Let CC be a control and communication law for S compatible with T. 

(i) The time complexity to achieve Twith CC from {xo,wo) G Hje/ "'^o ' ^ nie/ ^o*' '^^ 

TC(T, CC,xo, wq) = inf {£ \ T{x{tk).,w{tk)) = true , for all k> f\ , 

where 1 1— > {x{t),w{t)) is the evolution of {S,CC) from the initial condition {xq,wo). 

(a) The time complexity to achieve Twith CC is 

TC(r, CC) = sup { TC(r, CC,xo,wo) I (xo, wo) G J] 4'' x H } ■ 

(Hi) T/ie time complexity of T is 

TC(T) = inf{TC(75 CC) | CC is compatible with T} . • 
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Next, we define the notion of mean and totai communication complexities for a task. As usual, we 
assume that the network S has a communication edge map Ecmm and that the control and commu- 
nication law CC has language L and message-generation functions msgW, i G I. With these data we 
can discuss the communication cost of realizing one communication round. At time t G T from state 
(x, w) G Yliei ^ riie/ an element of L needs to be transmitted for each edge of the directed 
graph (/,Scmm\0(*,a;,'w)) defined by 

G Ecnim\<l){'t,x,w) if and only if G -Ecmm(a;) and msgW(t,3;W,'wW, j) / null. 

Next, we need a model for the cost of sending a message for each directed edge in -Ecmm\0- 

Definition 3.6 (One-round cost) (i) For / = {1, . . . , N}, a function Cj.^^ : 2^^^ — > M+ is a one- 
round cost function «/Cj,jj^(0) = 0, and 5*1 C 5*2 C / x / implies Cj.^^{Si) < Cj.j^^{S2)- 

(a) A one-round cost function Cj,^^ is additive if, for all Si, 52 C / x /, 5i n 5*2 = implies Cj.^j(5iU 
5*2) = Cj,„j(S'i) + Cj,^d(5'2). • 

We postpone our discussion about specific functions Cj-^j to the next subsection. Here we only empha- 
size that, for a given control and communication law CC with language L, the one-round cost depends 
on L; we therefore write it as C^^: 2^^^ — > M_|_. 

Definition 3.7 (Communication complexity) Let S be a robotic network and let T be a coor- 
dination task for S. Let CC be a control and communication law for S compatible with T, and let 
C^jj^: 2^^^ — > R+ be a one-round cost function. 

(i) Let {xq,wo) G riie/^o^ ^ riie/^o*'- mean communication complexity to achieve Twith 
CC from {xo,wo) and the total communication complexity to achieve Twith CC from (xo,wo) G 

Hie/ ^0 ' ^ Hie/ ^d*' '^^^> respectively, 

^ TC{r,CC,xo,wo)-l 

MCC{T,CC,Xo,Wo) = rj.r-,^^ rr T C^id ° ^cmm\0(*£, ^(t^)), 

l^[LL, 1,xq,wq) 

TC(r,CC,a;o,-(i;o)-l 

TCC(T,(:C,a;o,u;o) = ^ C^nd o ^cmm\0(i^, u;(t£)), 

where 1 1— > {x{t),w{t)) is the evolution of{S,CC) from the initial condition {xq,wq). (Here MCC 
is defined only for {xo,wo) with the property that T{xq,wo) = false. j 

(a) The mean communication complexity to achieve Twith CC and the total communication com- 
plexity to achieve Twith CC are, respectively, 

MCC{T,CC)=snp{MCC{T,CC,xo,wo) \ {xo,wo) G J] 4'' x H^d'')' 
TCC{T,CC) =snp{TCC{T,CC,xo,wo) \ (xcu'o) G H^o' ^ H^^^} • 

(Hi) The mean communication complexity of T and the total communication complexity of T are, 
respectively, 

MCC(T) = inf{MCC(T, CC) | CC is compatible with 7} , 

TCC(T) = inf{TCC(7; CC) \ CC is compatible with 7} . • 
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We conclude this subsection with some remarks. 

Remarks 3.8 (i) The total communication complexity is equal to the average transmission cost 
during the execution multiplied by the number of rounds required to carry out the desired task. 
That is, for {xo,wo) G Uiei ^0 x ^e/ 

TCC(T; CC,xo,wo) = MCC(7; CC,xo,wo) ■ TC(r, cc, xo, wo). 
In turn this implies that TCC(T,CC) < MCC(T,CC) • TC(T,CC). 

(ii) According to this notation, given a robotic network S and a control and communication law CC, 
the time complexity of achieving a task Twith CC is TC(T, CC) £ 0{f) (resp. TC(T, CC) E ^{f)), 
if there exist Nq £ N and k G such that TC(T,CC ,xq,wq) < kf{N) for all initial conditions 
{xq,wq) for each > Nq (resp. if TC{T,CC,xq,wo) > kf{N) for at least an initial condition 
{xo,wo) for each N > Nq). 

(iii) A different notion of communication complexity is defined in for a different robotic network 
model. Transcribed to the current setting, this notion of communication complexity of the 
execution of a control and communication law CC from initial conditions (xq, wq) would read as 

1 

cc(CC,xo,u;o)= lim -"^C^^^o Ec^^\^{ti,x(ti),w{t(>)) . (3) 

~pOO rx/ 

Note that this definition does not make reference to the completion of a task. We will later come 
back to this notion in Section f3. 41 • 

3.3 Communication costs in unidirectional and omnidirectional wireless channels 

In this subsection we discuss some modeling aspects of the one-round communication cost function 
/ X I D £" I— > Cj,jjjj(£') described in Definition l3.6l First, let us mention that the definition is motivated 
by the assumptions that (1) the cost of exchanging any message is bounded, and that (2) this cost is 
zero only for the null message. More specific detail about the communication cost depends necessarily 
on the type of communication service available between the agents. 

In unidirectional models of communication messages are sent in a point-to-point-wise fashion. 
Certain forms of communication, such as those based on the TCP-IP protocol, and certain technologies, 
such as wireless networks equipped with unidirectional antennas [22 . E3 | , fall into this category. On 
the other hand, in an omnidirectional model of communication (e.g., wireless networks equipped with 
omnidirectional antennas), a single transmission made by a node can be heard by several other nodes 
at the same time. This has the advantage that, by choosing a sufficiently large transmission power, a 
signal can reach all the neighboring nodes in a single time instant. 

Broadly speaking, it is very difficult to come up with an abstract model that captures adequately 
the cost of all possible communication technologies. For example, networking protocols for omnidi- 
rectional wireless networks rely on a many nested layers to handle, for example, media access, power 
control, congestion control, and routing. The presence of these layers and the non-trivial interactions 
between them make it difficult to assess communication costs of individual messages. Let us elaborate 
on this point in the following remark. 

Remark 3.9 (Omnidirectional wireless communication) The Minimum Power Broadcast (MPB) 
problem, the Medium Access Control (MAC) problem, and their relationship are subjects of vigorous 
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research in the wireless communications hterature, see for instance |2J| and references therein. Loosely 
speaking, the MPB problem consists of finding, for each agent i, the minimum broadcast radius i?W 
such that if agent i sends a message with communication radius ii'*', then its neighbors in a given 
graph E receive it. The MAC problem consists of determining a minimum number of broadcasting 
turns required for all agents to communicate their messages without interference. A schematic ap- 
proach to these problems is as follows: first, from the communication graph {I,E), one constructs the 
neighbor-induced graph (/, Ej\f) by 

{i,j) G Ej\f if and only if {i,j) £ E or (i, k), (j, k) £ E , for some k £ I 

In the new graph {I,Ej\f), the set of neighbors of the agent i is composed by its neighbors in the 
graph {I,E), together with the their respective neighbors. As a second step, one has to compute the 
chromatic number of the graph, i.e., the minimum number of colors xi^^Af) needed to color the agents 
in such a way that there are no two neighboring agents with the same color. (This is also referred 
to as the coloring- graph problem.) Theorem 5.2.4 in asserts that if a connected graph is neither 
complete, nor an odd cycle, then xi^Af) is less than or equal to the maximum valency of the graph. 
Once the chromatic number has been determined, broadcasting turns can be established according to 
an ordered sequence of the agents' colors. Although this approach is clearly inadequate, it provides 
some basic pointers with regards to communication costs. • 

Motivated by the difficulty of obtaining a detailed model, the rest of this paper relies on the 
following simplified models that capture some broad relevant aspects: 

(i) For a unidirectional communication model, E i— > Cj.^^{E) is proportional to the total number 
of non-null messages sent over the directed edges in E, that is, Cj.^^{E) = cq • cardinality (i?), 
where cq € M_|_ is the cost of sending a single message. This one-round cost function is additive. 
This number is trivially bounded by twice the number of edges of the complete graph, which 
is A'^(A'^ — 1). Therefore, for unidirectional models of communication, we have MCC(CC,7) G 
0{N^). 

(ii) For an omnidirectional communication model, E Cj.^^^{E) is proportional to the number of 
turns employed to complete a communication round without interference between the agents 
(see Remark This number is trivially upper bounded by N. Therefore, for omnidirectional 
models of communication, we have MCC(CC ,T) S 0{N). 

3.4 Invariance under rescheduling of control and communication laws 

In this section, we discuss the invariance properties of the notions of time and communication com- 
plexity under the rescheduling of a control and communication law. The idea behind rescheduling 
is to "spread" the execution of the law over time without affecting the trajectories described by the 
robotic agents of the network. There are at least two natural ways of doing this. One possible way 
consists of the network slowing down its motion, and letting some communication rounds pass without 
effectively interchanging any messages. Another possible way is to schedule the messages originally 
sent at a single time instant to be sent over multiple consecutive time instants, and adapt the motion 
of the network accordingly. Our objective is here is to formalize these ideas and to examine the effect 
that these processes have on the notions of complexity introduced earlier. For simplicity we consider 
the setting of static laws; similar results can be obtained for the general setting. 

Let S = (/, A, Ecmm) be a robotic network with driftless physical agents, that is, a robotic network 
where each physical agent is a driftless control system. Let CC = (No, I/, {ctlWjjg/, {msgWjjg/) be a 
static control and communication law. It is out intention to define a new control and communication 
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law by modifying CC; to do so we introduce some notation. Let s € N, with s < N, and let Vj = 
{lo, . . . , Is-i} be an s-partition of /, that is, Iq, . . . , /<j-i are disjoint and nonempty subsets of / and 

For i £ I, define the message-generation functions msg|'j : Nq x x I — > L by 



^^^m ^ /™sg['kiL^/sJ,a;,j), if i G 4 and A: = ^(mod s) , 

^ " I null, otherwise. 



According to this new message-generation function, only the agents with unique identifier in Ik will 
send messages at time t^, with £ G {k + as \ a £ Nq}. Equivalently, this can be stated as follows. Define 
the increasing function i^: Nq ^ Nq by F{£) = s{£ + 1) — 1. According to the message-generation 
functions specified by (@J), the messages originally sent at the time instant ti are now rescheduled to 
be sent at the time instants tp(^£^_s_^_i, . . . , tp^^^y Fig. 0] illustrates this idea. 




Figure 4: Under the rescheduling, the messages that are sent at the time instant ti under the control 
and communication law CC are rescheduled to be sent over the time instants . . . under 

the control and communication law CCi^g -p^y Accordingly, the evolution of the robotic network under 
the original law during the time interval [t^, t^+i] is now executed when all the corresponding messages 
have been transmitted, i.e., along the time interval 

For i G /, define the control functions ctl'^l : 1+ x x x ^ C/H by 

[i] J-^^^^l^^^-r^y^^ctlW (/i^(t),x,Xsmpid,y) , if t G and ^ = -l(mods) 



0, otherwise 



(5) 



where F -"^ : Nq — > Nq is the inverse of F, defined by F ^{i) = — 1, and for I = — l(mods), the 
function hi: [^^-^(f)) ^F-i(£)+i] is the time re-parameterization function defined by 

hi{t) = , t G [ti, tf+ij . 

Roughly speaking, the control law ctl|*j makes the agent i wait for the time intervals [t^, with 
i G {as — 1 I a G N}, to execute any motion. Accordingly, the evolution of the robotic network 
under the original law CC during the time interval t^+i] now takes place when all the corresponding 
messages have been transmitted, i.e., along the time interval [tF{i)-,tF{i)+i\- 
We gather the above construction in the following definition. 
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Definition 3.10 (Rescheduling of control and communication laws) Let S = {I,A,Ecmm) be 
a robotic network with driftless physical agents, and let CC = (Nq, {ctl'*l}jg/, {msgWjig/) be a static 
control and communication law. Let s G N, with s < N, and let Vi be an s-partition of I. The 
control and communication law CC (^g-p^) = (Nq, L, {ctl^*J p^^jjg/, {msg|*j -p^)}i6/) defined by equations @ 
and © is called a (s, P/)-rescheduling of CC. • 

Next, we examine the relation between the evolutions and the time and communication complexities 
of a control and communication law, and of those of its reschedulings. 

Proposition 3.11 With the same assumptions as in Definition \3.1IA let t i— > x{t) and t ^ x{t) 
denote the network evolutions starting from xq £ Hje/ ^ under CC and CC(-s -p^), respectively, and 
let T: Wiizj — > BooleSet be a coordination task for S. 

(i) For all k £ Nq, 



, for allt £ [tF(fc), iF(fc)+i] • 



(6) 



(a) For all Xq £ HiG/^o 

TC{CC(^s^-Pj),T,xo) = s • TC{CC,T,xo) . 

(Hi) //Cj.jjj is additive, then, for all xq £ flig/^o^ 

MCC(CC(,,p^),T,xo) = ^•MCC(CC,T,xo), 

and, therefore, the total communication cost of CC is invariant under rescheduling. • 

Proof: The relationships © are direct consequences of the definition of rescheduling. We leave 
the bookkeeping to the interested reader. By definition of TC(CC, TJ xq), we have that T{x{tk)) = true, 
for all k > TC(CC,T, xo), and T{x{t'YQ(^cc ,T,xo)-i)) — false. Let us rewrite these equalities in terms of 
the trajectories corresponding to the rescheduled control and communication law. From equation ©, 
one can write xW(tfc) = a;f*'(^F(A:)(iF(fe))) = ^^'^'(^^(fc))) for alH G I and k £ Nq. Therefore, we have 

r(x(t^(fc))) =T(x(tfc)) =true, for all F(A:) > F {TC (CC ,T, x^)) , 
T{x{tF{TC{CC,T,XQ)-i))) =T{x{tr^Qicc,T,xo)^i)) = false, 

where we have used the definition Q of the rescheduled message-generation function. Now, note that 
by equation one has 

= , for ah ^ e Nq and alH e J. 

Therefore, '?"(2;(tF(TC(CC,T,a;o)-i)+i)) = '^(^(*f(tc(CC,t,xo)))) we can rewrite the previous identities 
as 

T{i{tk)) = true , for all k > F(TC(CC, T, xq) - 1) + 1 , 
'T{x{tF{Tc{CC,T,xo)-i))) = false, 
which imply that 

TC(CC(,,p,),T,xo) = F(TC(CC,T,xo) - 1) + 1 = sTC(CC,T;xo) . 
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As for the mean communication complexity, additivity of Cj.^^ implies 

Crnd ° -E^cmm\0 ^cmm\0 

H h Crnd ° -^cmm\0 (tF(e),x{tF{e))) , 

where we have used F{i — 1) + 1 = F{i) — s + 1. Now, we compute 

TC(CC(,,p^),r,xo)-l F(TC(CC,r,xo)-l) 

e=o 1=0 

TC{CC,r,xo)-l F{i) TC{CC ,r,xo)-l 

= XI X] Cfnd ° -^cmm\0(*fc)^(*fc)) = X ^vnd° ^cram\ili{tliX{ti)) , 

i=Q k=F(i)-s+l 1=0 

which completes the proof of part (iii). ■ 

Remark 3.12 It is worth noting that the notion of communication complexity defined in is not 
invariant under rescheduling. Indeed, reasoning as before, one computes 

ks-l 



^™ h ^ ° ^cmm\0(*£, X{ti)) = _ lim -= ^ ^rnd ° Ecmm\d{ti, x{t(,)) 

^ i=0 ''^+°° ^■^ ~ 1 e=o 

^ k-1 m 

= _lim J "X X C^nd° E^mm\<i){ti,x{te)) 

fc— >+co KS — i - 

i=0 i=F{i-l)+l 

= - 1™ T 7 X ^md ° -^cmm\0 • 

fc^+oo ks — I - 



Therefore, cc(CC(s -p^), xq) = ^ cc{CC,xo). This means that, by performing a rescheduling of the control 
and communication law, one can indeed lower the measure of communication complexity cc, although 
the trajectory described by the robotic network will continue to be the same. • 

4 Motion coordination algorithms and their time complexity 

In this section we provide examples of motion coordination algorithms for robotic networks performing 
a variety of distributed tasks. For each algorithm and task, we present some results on the correspond- 
ing time and communication complexity. 

4.1 Agreement on direction of motion and equidistance 

From Examples 12.71 l^iTH and recall the definition of uniform network 5gi ^-disk of locally-connected 
first-order agents in S^, the agree-and-pursue control and communication law CCagr-pursuit > and the two 
coordination tasks Tdrctn and T^.cqdstnc- 

Theorem 4.1 For kp^op £]0, \[, r g]0, 27r], a = Nr — 27r and e g]0, 1[, the network Sgi f.^isk, the law 
CCagr-pursuit, and the tasks T^rctn and 7;_eqdstnc together satisfy: 

(i) the upper bound TC(7^rctn,CCagr-pursuit) G 0{Nr^^) and the lower bound 
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\n{N) if a < 0; 



(ii) if a > 0, then the upper bound TC(7^_eqdstnc5 CCagr-pursuit) G 0{N^ log{Ne ^) + A^loga ^) and 

the lower bound TC(7^_cq(jstnc;CCagr-pursuit 

) G n{N^log{e-^)). Ifa<0, then CC^^, 

-pursuit does 

not achieve 7^-cqdstnc in general. 

Proof: Let us start by proving fact (i). Without loss of generality, assume drctii[^l(0) = c, and 
that Tdrctn is false at time 0. Therefore, at least one agent is moving counterclockwise at time 0, and 
we can define k = max{i G / | drctii[*l(0) = cc}. Define tk = inf({^ G No | drctnW(^) = c} U {+oo}). 
In what follows we provide an upper bound on tk- 
For i < tk, define 

j{£) = argmin{distc(0W(^),^[''l(^)) | priori*^ = N, i e I}. 

In other words, for all instants of time when agent k is moving counterclockwise, the agent j{l) has 
prior equal to A^, is moving clockwise, and is the agent closest to agent k with these two properties. 
Clearly, 

27r > distc(0[^l(O),^'^ko)) = distc(0[-''(°^ko),0t''l(O)). 
Additionally, for ^ < — 1, we claim that 

distc(0[^'(^)]W,e[^'(^+i)l(^+ 1)) > kpropr. 

This happens because either (1) there is no agent clockwise-ahead of 6^^^^^\£) within clockwise distance 
r and, therefore, the claim is obvious, or (2) there are such agents. In case (2), let m denote the agent 
whose clockwise distance to agent is maximal within the set of agents with clockwise distance r 
from 6'[jW1(£). Then, 

distc{9^^^^'>\£),e^^'^^+^^^{£ + l)) = distc(^[^'(^)l(^),^^™k^ + 1)) 

= distc(^[^'(^)lW,^^™k^)) +distc(e['"l(^),^[™J(^ + 1)) 

> distc(^[^'(')lW,^'™k^)) + kprop{r - dist,(0[^'W](^),^H(^))) 

= kp,opr + (1 - fcprop) dist,(0bW] (^), (^)) > k^^^^r, 

where the first inequality follows from the fact that at time I there can be no agent whose clockwise 
distance to agent m is less than (r - distc(6'y(^)l(£), (£))). 

In summary, either agent k changes direction of motion or at each instant of time its distance to 
the closest agent with prior equal to A^ decreases by a constant feprop^. This shows that 

tk G 0{r~'). 

Because at least one more agent moves in the clockwise direction after 0{r~^) time, it follows that all 
agents will move clockwise after 0{Nr~^) time. This completes the proof of the upper bound in (i). 
Let us prove the lower bound in (i). Note that the number of agents that can fit into the circle 
spaced a distance r apart one from each other, is at most [^J . Consider an initial configuration where 
drctnW(O) = cc for i G {1, . . . , A^ - 1}, drctn[^l(0) = c, and 

(i) for a > 0, \_^\ agents (including the agent A^) lie at a distance r one from each other in S^, and 
the remaining A^ — [^J agents are located within a clockwise distance | of the agent i* lying in 
the position in symmetric to the location of agent A^; 

(ii) for a < 0, all agents lie at a distance r in in counterclockwise order incrementally according 
to their unique identifier (note that they might not cover the whole circle). 



20 



Note that the displacement of each agent is upper bounded by kp^-opr < |. When a > 0, the number 
of time steps that takes agent i* to receive the message with priority N is lower bounded by ^ . 
Therefore, TC(7^rctn, CCagr-pursaiiiJ ^ 0.{r~^). When a < 0, the analysis of the network evolution is 



parallel to the one discussed in |27l . Chapter 1] for the leader election algorithm in static networks with 
the ring topology. The number of time steps that takes agent 1 to receive the message with priority 
is given by iV - 1. Therefore, TC(Tdrctn,CCagr-pursuit) £ ^{N). 

To prove fact (ii), we assume that T^rctn has been achieved (so that all agents are moving clockwise), 
and we first prove a fact regarding connectivity. At time £ G No, define 

H{i) = {x G I mindistc(a;,^[^](£)) + mindistcc(a;, > r}. 

i&i j&i 

In other words, any point in H{t) is at least a distance r, clockwise or counterclockwise, from an 
agent. Therefore, H{€) does not contain any point between two agents separated by a distance less 
than r, and each connected component of H{1) has length at least r. Let nni^) be the number of 
connected components of H(i), if H(£) is empty, then we take the convention that nni^) = 0. Clearly, 
nni^) ^ N ■ We claim that, if nni^) > 0, then t i-^ + t) is non-increasing. Let d{£) < r be 

the distance between any two consecutive agents at time i. Because both agents move in the same 
direction, a simple calculation shows that 

d{£ + 1) < d{£) + kprop{r - d{£)) = (1 - kpT.op)d{£) + kp^o^r < (1 - A:prop)r + /cpropr- = r. 

This means that the two agents remain within distance r and, therefore connected, at the following 
time instant. Because the number of connected components of Er{d^^\ . . . , 6^^^) does not increase, it 
follows that the number of connected components of H cannot increase. 

Next we claim that, if n}{{£) > 0, then there exists t > £ such that nnit) < n}{{£)- By contradiction, 
assume nH{£) = nnit) for all t > £. Without loss of generality, let {!,..., m} be a set of agents with 
the properties that distcc {9^^ {£) , 6^'+^^ {£)) < r, for i e {l,...,m}, that 6'W(^) and 6'H(£) belong to 
the boundary of H{£), and that there is no other set with the same properties and more agents. One 
can show that, for t > £, 

0W(^t + l) = em{t)-kp,opr, 

0[^^t + 1) = e^Ht) - fepropdist,(0W(t),0[^-i](O), ie{2,...,m}. 

If we define d{t) = {distcc{0^^\t),9^^^{t)), . . . ,distcciO^"'~^\t),0^"'^{t))) G M™"\ then one can show 
that 

kr) 



d{t + 1) = Tridm_i(/cpi.op, 1 - kprop, 0) d{t) 



^prop 





where the linear map (a, 6, c) ^ Tridm-i(a, ^, c) G ]R('^~^)^(™'~^) is defined in Appendix^ This is 
a discrete-time linear time-invariant dynamical system with unique equilibrium point r(l, . . . , 1). By 
Theorem IB.3f ii) in Appendix El for ry g]0, 1[, the solution t ^ d(t) to this system reaches a ball of 
radius r] centered at the equilibrium point in time 0{mlogm + log?7~^). (Here we used the fact that 
the initial condition of this system is bounded.) In turn, this implies that t i— > ^^^idi^t) is larger 
than (m — l)(r — r]) in time 0(mlogm + \ogi]~^). 

We are now ready to find the contradiction and show that nnit) cannot remain equal to nH{£) for 
all time t. After time 0(m log m + log ry"-"^) = 0(A^log + logry"-*^), we have: 

2tt > nH{£)r + ^ {r - 'n){mj - 1) = nH{£)r + {N - nH{£)){r - rj) = nH{£)v + N{r - rj). 

3=1 
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Here mi, . . . , mnjj(i) are the number of agents in each isolated group, and each connected component 
of H(i) has length at least r. Now, take ij 



Nr-2n 
N 



and the contradiction follows from 



27r > nuitfr] + Nr - Nr] = nuitfr] + Nr + 2tt - Nr = nuitjT] + 27r. 

In summary this shows that, in time 0(A^log + log 77^"*^) = 0(A^log + loga~^), the number 
of connected components of H will decrease by one. Therefore, in time 0(A^^logA^ + A^loga"^) 
the set H will become empty. At that time, the resulting network will obey the discrete-time linear 
time-invariant dynamical system: 

d{t + 1) = CircAr(A;pi.op, 1 - Vop, 0) d{t). 

Here d{t) = ( distcc(^W(t), ^'^l (*))>•••, distcc(^[^l(*)) ^'^^^'(0)) e with the convention ^t^+i] = 
6l[^l. By Theorem ESJiii) in 

Appendix El the solution t ^ d{t) reaches the desired configuration in 
time 0{N'^ loge~^) with an error whose 2-norm, and therefore, its 00-norm is of order e. In summary, 
the desired configuration is achieved in time 0{N'^ log(A^e~^) -|- log a~"^). 

For the lower bound, consider an initial configuration with the properties that (i) agents are 
counterclockwise-ordered according to their unique identifier, (ii) the set H is empty, and (iii) the 
inter-agent distances d{0) = ( distcc(^W (0), ^^^1 (0)), ■ ■ ■ , distcc(e[^l(0), (0))) are given by 



d(0) 



2n' 

N 



2n 

.N. 



where vjy is the eigenvector of CircAr(A;prop, 1 — fcpropjO) corresponding to the eigenvalue 1 — /cprop + 
^propcos (^) — A;prop\/— Tsin (^) (see Appendix 0, and > is chosen sufficiently small so that 
d{0) S M.^. By Theorem IB . Sr iii) in Appendix^ the solution 1 1— > d{t) reaches the desired configuration 
in time G(A^^loge~^) with an error whose 2-norm, and therefore, its 00-norm is of order e. This 
concludes the result. ■ 



4.2 Rendezvous without connectivity constraint 

From Examples 12.41 HTT^ and recall the definition of uniform network S^d ,,_disk °f locally-connected 
first-order agents in M"', the move-toward-average control and communication law CCavrg, and the 
coordination task T^ndzvs- 

Theorem 4.2 For d = 1, the network S^d r_disk' ^'^^ CCavrg? and the task T^ndzvs together satisfy 

TC(7;.„dzvs,CCavrg) G 0{N^) and TC(T,Mzv1, CCavrg) G ^{N). 

Proof: One can easily prove that, along the evolution of the network, the ordering of the agents 
is preserved, i.e., if < x^\l), then x^^\£+l) < ^[■''(^-l- 1). However, links between agents are not 

necessarily preserved (see e.g. Figure^. Indeed, connected components may split along the evolution. 
However, mergings are not possible. Consider two contiguous connected components Ci and C2, with 
Ci to the left of C2. By definition, the rightmost agent of Ci and the leftmost agent of C2 are at 
a distance strictly bigger than r. Now, by executing the algorithm, they can only but increase that 
distance, since the rightmost agent of Ci will move to the left, and the leftmost agent of C2 will move 
to the right. Therefore, connected components do not merge. 

Consider first the case of an initial configuration of the network for which the communication graph 
remains connected throughout the evolution. Without loss of generality, assume that the agents are 
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ordered from left to right according to their identifier, that is, a;[^](0) = {xq)i < ■ ■ ■ < x[^](0) = (xo)Ar- 
Let a £ {3, ... , N} have the property that agents {2, . . . , a — 1} are neighbors of agent 1, and agent 
a is not. (If instead ah agents are within an interval of length r, then rendezvous is achieved in 1 
time instant, and the statement in theorem is easily seen to be true.) Note that we can assume that 
agents {2, . . . , a — 1} are also neighbors of agent a. If this is not the case, then those agents that are 
neighbors of agent 1 and not of agent a, rendezvous with agent 1 at the next time instant. At the 
time instant i = 1, the new updated positions satisfy 

^ Q— 1 ^ a 

xW (1) = (0), ^'^1 (1) e [- E (0), *1 , 7 G {2, . . . , a - 1}, 

fc=i fc=i 

where * denotes certain unimportant point. 
Now, we show that 

xW (a - 1) - (0) > — -. (7) 

aya — 1) 

Let us first show the inequality for a = 3. Note that the fact that the communication graph remains 
connected implies that agent 2 is still a neighbor of agent 1 at the time instant i = 1. Therefore 
2;[^1(2) > + and from here we deduce 

xW(2) - xW(0) > - xW(0)) 

> ^(^(xW(0) +xPl(0) +x[3l(0)) -xW(0)) > i(x[3l(0) -xW(0)) > ^. 

Zi o DO 

Let us now proceed by induction. Assume that inequality Q is valid for a — 1, and let us prove it for 
a. Consider first the possibility when at the time instant i = I, the agent a — 1 is still a neighbor of 
agent 1. In this case, xW(2) > Ylk=i x^''^^), and from here we deduce 

xW(2) - ^™(0) > - xW(0)) > y xW(0) - xW(0)) 



k=l 



> . ^_ (x["l (0) - xt^l (0) ) > 



a{a — 1) V / a{a — 1) ' 

which in particular implies 0. Consider then the case when agent a — 1 is not a neighbor of agent 
1 at the time instant i = 1. Let (3 < a such that agent /3 — 1 is a neighbor of agent 1 at i = 1, but 
agent /3 is not. Since /3 < a, we have by induction x^^\(3) — ^[^^(l) > j^^^^zj^- From here, we deduce 
thatxW(a-l)-xW(0)>^. 

Inequality implies that, at most ina — 1<A^ — 1 time instants, the leftmost agent traverses a 
distance greater than or equal to (provided that at each step there exists at least another agent 

which is not its neighbor). Since diam(xo, /) < {N — l)r, we deduce that in A^(A^ — 1)'^ time instants 
there cannot be any agent which is not a neighbor of the agent 1. Hence, all agents rendezvous at the 
next time instant. Consequently, 

TC(rrndzvs,CCavrg,Xo) < N{N - if + 1. 

Finally, for a general initial configuration xq, because there are a finite number of agents, only a finite 
number of splittings (at most — 1) of the connected components of the communication graph can 
take place along the evolution. Therefore, we conclude TC(7^ndzvs5 CCavrg) = 0{N^). 
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Let us now prove the lower bound. Consider an initial configuration xq G M where all agents 
are positioned in increasing order according to their identity, and exactly at a distance r apart, say 
(xo)j+i — {xo)i = r, i £ {1, . . . , N — 1}. Assume for simplicity that N is odd - when N is even, one can 
reason in an analogous way. Because of the symmetry of the initial condition, in the first time step, 
only agents 1 and move. All the remaining agents remain in their position because it coincides with 
the average of its neighbors' position and its own. At the second time step, only agents 1, 2, — 1 
and A^ move, and the others remain still because of the symmetry. Applying this idea iteratively, one 
deduces the time step when agents ^^^-^ and ^^-^ move for the first time is lower bounded by ^^^y^- 
Since both agents have still at least a neighbor (agent ^-^), the task T^ndzvs has not been achieved 
yet at this time step. Therefore, TC(7^ndzvs)CCavrg, a^o) > ^^^y^, and the result follows. ■ 

4.3 Rendezvous with connectivity constraint 

In this section we shall consider both networks S^d ^.^jg^ ai^d S^d presented in Examples 12.41 
andEH 

Circumcenter control and communication law 

Here we define the circumcenter control and communication law CCcrcmcntr for both networks S^d ^.^isk 
and S^d r-i£L- This is a uniform, static, time-independent law originally introduced by Q] and later 
studied in |5|, |25| . Loosely speaking, the evolution of the network under the circumcenter control and 
communication law can be described as follows: 

[Informal description] Communication rounds take place at each natural instant of time. 
At each communication round each agent performs the following tasks: (i) it transmits 
its position and receives its neighbors' positions; (ii) it computes the circumcenter of the 
point set comprised of its neighbors and of itself, and (iii) it moves toward this circumcenter 
while maintaining connectivity with its neighbors. 

Let us present this description in more formal terms. We set T = No, L = M^, and msgW = msgg^^j, 
i £ I. In order to define the control function, we need to introduce some preliminary constructions. 
First, connectivity is maintained by restricting the allowable motion of each agent in the following 
appropriate manner. If agents i and j are neighbors at time £ G Nq, then we require their subsequent 
positions to belong to B{- — — |). If an agent i has its neighbors at locations {qi, . . . ,qi} at 
time £, then its constraint set T^^iii^g-j ({<?!) ■ ■ ■ ,Qi}) is 



Second, in order to maximize the displacement toward the circumcenter of the point set comprised 
of its neighbors and of itself, each agent solves a convex optimization problem that can be stated in 
general as follows. For qq and qi in M"^, and for a convex closed set Q CM.'^ with go S Q, let X{qo, qi, Q) 
denote the solution to the strictly convex problem: 



Under the stated assumptions the solution exists and is unique. Third, note that since the agents 
operate with the standard message-generation function, it is clear that the natural projection ii^d 




maximize A 



subject to A < 1, (1 — A)g'o + Agi G Q. 
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maps the messages received at time £ G Nq by the agent i G / onto the positions of its neighbors. 

We are now ready to define the last constitutive element of CCcrcmcntr- Define the control function 

ctl : R'^ X R'^ X ^ R'^ by 

ctl(j;, 2;smpld, y) = K - (Circum(7riRd(?/) U {xsmpld}) - a^smpld) , (8) 

where A* = A(xsmpid, {Civcum{Tr^d{y)U{xsmpid}),'Dx,^pia,riTrRd{y))). Evolving under this control law, it 
is clear that, at time [ij + 1, each agent i reaches the point (1 — A,,, )xW ( [tj ) + A* Circum(7rigd (y W ( [fj )) U 

{xW(LtJ)}). 

Next, we consider the network tSr-oo-disk 

of locally-oo-connected first-order agents in M'^, see Ex- 
ample 12.61 For this network we define the parallel circumcenter law, CCpU-crcmcntn by designing d 
decoupled circumcenter laws running in parallel on each coordinate axis of M'^. As before, this law is 
uniform, static and time- independent. As before, we set T = No, L = R'^, and msg[*l — msgstfj, i G I. 
We define the control function ctl : M'^ x M*^ x ^ M'^ by 

ctl(3;,2;smpid,y) = (Circum(ri(A4)) - ri(a;smpid), . . . , CiT:cum{Td{M)) - Tdix^nipld)) , (9) 

where M. = ■7r^d{y) U {a^smpid}, and ri, . . . ,rrf: M'^ — > M denote the canonical projections of onto R. 
See Fig. 15] for an illustration of this law in . 










T '■■ 

xl'i^ ' ~ ' 





Figure 5: Parallel circumcenter control and communication law in R^. The target point for the agent 
i is plotted in light gray and has coordinates (Circum(ri(7WW)), Circum(r2(A^ '*'))). 

Asymptotic behavior and complexity analysis 

The following theorem summarizes the results known in the literature about the asymptotic properties 
of the circumcenter control and communication law. 

Theorem 4.3 (Correctness of the circumcenter law) For d G N, r G M_|_ and e G M+, the fol- 
lowing statements hold: 

(i) on the network S^d j._disk' ^'^^ CCcrcmcntr achieves the exact rendezvous task T^ndzvs,' 

(a) on the network S^d ^-lD' ^^'^ CCcrcmcntr achieves the e-rendezvous task T^-mdzvs/ 

(Hi) on the network S^d ,,_oo-disk' ^'^^ CCpii_crcmcntr achieves the exact rendezvous task T^ndzvs/ 
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(iv) the evolutions of {Sj^d j._^^^y.,CC crcmcntv), of {S^d^^_-^Y),CC cvcmcntr), and of {S^d CC pW.cvcmcntr) 

have the property that, if two agents belong to the same connected component of the communica- 
tion graph at some time £ G Nq, then they continue to belong to the same connected component 
of the communication graph for all subsequent times k > i. • 

Proof: The results on 5r-disk appeared originally in 0]. The proof for the results on Sr-LD is 
provided in We postpone the proof for Sr-oo-disk to the proof of Theorem 14.41 below. ■ 
Next we analyze the time complexity of CCcrcmcntr • We provide complete results only for the case 
d = 1. As we see next, the complexity properties of CCcrcmcntr differ dramatically when applied to the 
two robotic networks with different communication graphs. 

Theorem 4.4 (Time complexity of circumcenter law) For r G M+ and e g]0, 1[, the following 
statements hold: 

(i) for d=l, on the network cSM,r-disk, TC(7;ndzvs, CCcrcmcntr) G Q{N); 

(a) ford= 1, on the network Sm^r-LD, TC(T(re)_i.ndzvs, CCcrcmcntr) G @{N'^log{Ne~^)); 

(Hi) for d£N, on the network S^d TC(Trndzvs, CCpii_crcmcntr) G @{N). • 

Proof: Let xq G M^. Throughout the proof, neighboring relationships are understood with 
respect to the r-disk graph. First of all, let us show that, for n = 1, the connectivity constraints 
on each agent i £ I imposed by the constraint set ^('/rK(y))) are superfluous, i.e., the solution 
of the convex optimization problem is A* = 1 (cf. equation ©). To see this, assume that agents 
i and j are neighbors at time instant £, define A^W as TTj^d (£)) U {xW(^)}, and let us show that 
Circum(A^W) belongs to ;g( ^''' W+^'^'ffl ^ Without loss of generality, let x^^{e) < x^^^i). Let x^^{e), 
x^^{£) denote the positions of the leftmost and rightmost agents among the neighbors of agent i. Note 
that xH(^) < x^j\i) < x^^{e) and Circum(AlW) = ^{x^^{i) + x^^{e)). Then, 

I Circum(M[^l) - ^(xW(£) + x^\£))\ = ^\x^^{i) - x^^{e) + x^^{£) - x^\£)\ 

< ^ma^{\x^^{£) - x^'\£)\,\x^^{i) - x^^\£)\} < ^, 

as claimed. Therefore, we have that x[*l(£+l) = Circum(AlW). Likewise, one can deduce Circum(AlW) < 
Circum(A^[-']), and therefore, the order of the agents is preserved. 

Consider first the case when ^^r-disk(3^o) is connected. Without loss of generality, assume that the 
agents are ordered from left to right according to their identifier, that is, ^[^^(O) = (xo)i < ••• < 
a;[-^](0) = {xo)n- Let a G {3, . . . , A^} have the property that agents {2, ... ,a — 1} are neighbors of 
agent 1, and agent a is not. (If instead all agents are within an interval of length r, then rendezvous 
is achieved in 1 time instant, and the statement in theorem is easily seen to be true.) See Fig.Elfor an 
illustration of these definitions. Note that we can assume that agents {2, . . . , a — 1} are also neighbors 
of agent a. If this is not the case, then those agents that are neighbors of agent 1 and not of agent 
a, rendezvous with agent 1 at the next time instant. At the time instant £ = 1, the new updated 
positions satisfy 

J](l),^'"(0)+_^=--'-"(0), ,H(i) 



x W (0) + (0) xt^l (0) + (0) + r 



, 7e {2, 



, a 



!}• 
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Figure 6: Definition of a G {3, . . . , A^} for an initial network configuration. 
These equalities imply that — x[^](0) = ^(x[°"^](0) — ^[-^^(O)) < ^r. Analogously, we deduce 



3;[^1(2) — x[^l(l) < ^r, and therefore 



xW(2) -xW(0) < r. (10) 

On the other hand, from x^^\2) G [^(a;l^](l) + x^""^' (1)) , *] (where the symbol * represents a certain 
unimportant point in M), we deduce that 

a;W (2) - (0) > ^ (x^ (1) + xl"-^! (1)) - x™ (0) > ^ (xl^^^] (1) - x™ (0)) 

> i(^M±^ _ = 1(,H,„, _,(.!(„)) > (u, 

Inequalities (|lflj) and (|1H) mean that, after at most two time instants, agent 1 has traveled an amount 
larger than r/4. In turn this implies that 

diam(xo,/) 4diam(xo, J) 

^ !b Lj^^ij.n<j2vs 5 '-''-'crcmcntr ) '';o J _i ^ 

If -Er-disk(3;o) is not connected, note that along the network evolution, the connected components of 
the r-disk graph do not change. Therefore, using the previous characterization on the amount traveled 
by the leftmost agent of each connected component in at most two time instants, we deduce that 

1 4 

max diam(xo,C) < TC(7;ndzvs,CCcrcmcntr,a;o) < - max diam(xo,C). 

Note that the connectedness of each C £ C£;^_^j^^(xo) implies that diam(xo,C) < {N — l)r, and 
therefore TC(7^ndzvs)CCcrcmcntr) £ 0{N). Moreover, for xq G such that (xo)i+i — (xo)j = r, 
z G {1, . . . , — 1}, we have diam(xo, /) = {N — l)r, and therefore TC(7^ndzvs5 CCcrcmcntr, 2:0) > — 1- 
This concludes the proof of fact (i): 

TC('7^n(j2vs, CCcrcmcntr) ^ 0(iV) . 

Next we prove fact (ii). In the r-limited Delaunay graph, two agents on the line that are at most at 
a distance r from each other are neighbors if and only if there are no other agents between them. Also, 
note that the r-limited Delaunay graph and the r-disk graph have the same connected components 
(cf. 0). Using an argument similar to the one above, one can show that the connectivity constraints 
imposed by the constraint sets set ^^a;I'l([tJ),r(^iR(y))) are again superfluous. 
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Consider first the case when -Br-Lola^o) is connected. Note that this is equivalent to £'r-disk(2^o) 
being connected. Without loss of generality, assume that the agents are ordered from left to right 
according to their identifier, that is, x M(0) = (2:0)1 < ••• < x[^](0) = (xo)jv. The evolution of the 
network under CCcrcmcntr can then be described as the discrete-time dynamical system 



1 



4^-i](£+ 1) = l(^[^-2](^) ^^[JV](^)) x^^\£+l) = + 



Note that this evolution respects the ordering of the agents. Equivalently, we can write x{i + 1) 
Ax{£), where A is the N x N matrix given by 



A 
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Now, note that A = ATrid^(^,0) as defined in Appendix IHl Theorem IB.4r i) implies that, for 
^ave ~ j^^"^ "^{^^ 1 ^6 have that lim£_>_|_oQ x{i) = Xavel, and that the maximum time required for 



\x{e) 



Xr. 



< r]\\x{0) — Xa 



(over all initial conditions x(0) £ M^) is 0(A^^log77 (As an 



aside, this also implies that the agents rendezvous at the location given by the average of their initial 
positions. In other words, we can forecast the asymptotic rendezvous position for this case, as opposed 
to the case with the r-disk communication graph.) 

Next, let us convert the contraction inequality on 2-norms into an appropriate inequality on 00- 
norms. Note that diam(a;o,/) < {N — l)r because Er-LDi^o) is connected. Therefore 



\x{0) 



X. 



avc -^00 



max 1x^*^(0) 



[1] JN] 







X, 







< (iV - l)r. 



For £ of order iV^logry ^, we use this bound on ||x(0) 
||f II2 < \/iV||f lloo for all V G M^, to obtain: 



Cavcl||oo and the basic inequalities ||f ||oo < 



\x{i) 



f-avel II 00 



< ||x(£) - Xavellb < vWxiO) 



< riVN\\x{0) - Xavcllloo < vVN{N - l)r. 



This means that (re)-rendezvous is achieved for 7]^/N{N — l)r = re, that is, in time 0(A^^log?/ ^) = 
0{N^ log(iVe-i)). 

Next, we show the lower bound. From equation 1)15(1 in the proof of Theorem IB. 31 we recall the 
unit-length vector v^v-i G This vector is an eigenvector of Tridjv-i(^, 0, |) corresponding 



to the largest singular value cos(-^). For ^ 



IOV2 



'jV^/^, we then define the initial condition xq 



VAT-l 



G R . One can show that (xo)i < (xo)j+i for i G {1, . . . , A^— 1}, that (xq) 



0, and that 



max{(xo)j+i — {xo)i | i G {1, . . . , A — 1}} < r. By Lemma FB . 51 and because 



\w\ 



< \\w\\2 < ^/NWwI 



for all w G 



, we compute 

^^5/2 I 



10^2 







00 ~ 10\/2 







2 ~ 10^2 



VAr_i||2 



rA 
IOV2' 
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The trajectory x{tj = (cos(;^))^3;o therefore satisfies 

\W)\U = (cos (^)) Volloo > (^))'- 

Therefore, ||x(^)||oo is larger than ^re so long as -j^^-j^N {cos{j^)Y > ^e, that is, so long as 

^ ^ log(g^i iV) - log(5V2) 
-log(cos(f)) 

The rest of the proof is the same as in Theorem IB.Sf i) for the lower bound result. 

If Ej._ijq{xq) is not connected, note that along the network evolution, the connected components 
do not change. Therefore, the previous reasoning can be applied to each connected component. Since 
the number of agents in each connected component is strictly less that N , the time complexity can 
only but improve. Therefore, we conclude that 

TC(T,ndzvs,CCcrcmcntr) G Q{N^ \og{N e~^)) . 

This completes the proof of fact (ii). 

Finally, we prove the statements regarding S^d r_oo-disk ^'^^'^ CCpu.crcmcntr in fact (iii) and in the 
previous Theorem 14.31 By definition, agents i and j are neighbors at time £ G No if and only if 
||a;W(^) — a;[-'l(^)||oo < r, which is equivalent to 

|rfc(xWW)-Tfc(xb']W)l <r, A:G{l,...,d}. 

Recall from the proof of fact (i) that the connectivity constraints of CCcrcmcntr on each agent are 
trivially satisfied in the 1-dimensional case. This fact has the following important consequence: from 
the expression for the control function in CCpU-crcmcntrj we deduce that the evolution under CCpu_crcmcntr 
of the robotic network S^d ^.oo.disk (in d dimensions) can be alternatively described as the evolution 
under CCcrcmcntr of d robotic networks 5R,r-disk in IR- The correctness and the time complexity results 
now follows from the analysis of CCcrcmcntr at d = 1. ■ 
Finally we proceed to characterize the mean communication complexity of the circumcenter control 
and communication law. We consider the case of a unidirectional communication model with one-round 
cost function depending linearly on the cardinality of the communication graph. 

Theorem 4.5 (Mean communication complexity of circumcenter law for unidirectional communicatior 

For r G and e g]0, 1[, the following statements hold: 

(i) for d=l, on the network S^d^^_^^^y., MCC(7;ndzvs, CCcrcmcntr) G @{N'^); 

(ii) for d= 1, on the network cSjjd MCC(7;ndzvs, CCcrcmcntr) G &iN); 

(iii) for d£N, on the network S^d ^-oo-disk' MCC(7;ndzvs, CCpu 

-crcmcntr ) G e{N^). • 

Proof: We start by proving fact (i). Let xq G (M'^)^ be such that all (xo)i, i & I, belong to a 
closed ball of radius \/^r/A. In such a case, one can deduce that (i) -Er-disk(a;o) is the complete graph, 
and therefore all agents compute the same goal point Circum, and (ii) this circumcenter belongs to the 
constrained set of each agent i £ I. As a consequence, all mobile agents rendezvous at the same location 
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Circum in one time instant. For such initial conditions, we have TC(7^ndzvs) CCcrcmcntr, 2:0) = 1, and 
therefore 

^ TC(7^ndzvs!^^crcmcntri^o) 1 

MCC(7;ndzvs, CCcrcmcntr, a;o) = 7f^=w7f JTp T X] ° -E'cmm\0(^, 

rndzvs ) '-' '-- crcmcntr ) ^0 J 



1=0 



N{N - 1) 



This proves that MCC(7^ndzvs;CCcrcmcntr) G 0(A^^), that is, fact (i). Additionally, fact (iii) is proved 
by analyzing the parallel circumcenter law as d decoupled versions of the circumcenter law. Next we 
show fact (ii). Given that the r-limited Delaunay graph is a subgraph of the Delaunay graph 0], 
and the number of edges of the latter is bounded by 3A^ — 6 (see, for instance, we deduce that 

MCC(7^ndzvs) CCcrcmcntr) S 0{N). On the other hand, for any initial configuration xq where all the 
agents are aligned and non-coincident, there are — 1 neighboring relationships which are preserved 
throughout the network evolution. In this case, MCC(7^ndzvs)CCcrcmcntr; a^o) = 2(A^— 1), and the result 
in fact (ii) follows. ■ 

Remark 4.6 Theorems 14.41 and 14.51 induce lower bounds on the time and mean communication com- 
plexity of the circumcenter law for the higher-dimensional case. Indeed, as a consequence of these 
results, we have 

(i) for d G N, on the network cSR^^.disk, TC(7;ndzvs, CCcrcmcntr) G ^{N) and MCC(7;ndzvs, CCcrcmcntr) e 

(ii) for d G N, on the network cSR,r--LD, TC(T(rj)_i.ndzvs, CCcrcmcntr) G il(iV^ log(iVe~^)) and MC C (T^ndzvs, CC 

n{N). 

We have performed extensive numerical simulations for the case d = 2 and the network S^d ^_disk- 
We have ran the algorithm starting from generic initial configurations (where, in particular, agents' 
positions are not aligned) contained in a bounded region of M?. We have consistently obtained that 
the time complexity to achieve T^-ndzvs with CCcrcmcntr starting from these initial configurations is 
independent of the number of agents. This leads us to conjecture that, in fact, initial configurations 
where all agents are aligned (i.e., the 1-dimensional case) give rise to the worst possible performance 
of the algorithm. In more formal terms, we conjecture that, for d > 2, TC(7^ndzvs) CCcrcmcntr) = S(-A). 



4.4 Deployment 

In this section we consider the uniform robotic network S^d presented in Example 12.51 with pa- 
rameter r G ]R+. We assume we are given a convex simple polytope Q C M'^, with an integrable density 
function : Q — > M_|_ . We assume that the initial positions of the agents belong to Q and we intend to 
design a control law that keeps them in Q for subsequent times. To achieve the e-r-deployment task 
discussed in Example 13.41 we define the centroid control and communication law CCcrcmcntr- This is a 
uniform, static, time-independent law studied in 0, 0| ■ Loosely speaking, the evolution of the network 
under the centroid control and communication law can be described as follows: 

[Informal description] Communication rounds take place at each natural instant of time. 
At each communication round each agent performs the following tasks: (i) it transmits its 
position and receives its neighbors' positions; (ii) it computes the centroid of an appropriate 
region (the region is the intersection between the agent's Voronoi cell and a closed ball 
centered at its position and of radius |), and (iii) it moves toward this centroid. 
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Let us present this description in more formal terms. We set T = No, L = W^, and msg^ = msgg^^j, 
i e I. We define the control function ctl : M'^ x M'^ x ^ R'^ by 

ctl(x,a;smpid,y) = Centroid(A') - Xs^pid , 

where X = Q S(xsmpid, §) n ( npe^^(j^) i^x,,„pM,p) and i^x^n^pid.P is the half-space {g e M'^ | \\q - 
f^smpidlb < Ik — Plb}- One can show that is a positively-invariant set for this control law. 

The following theorem on the centroid control and communication law summarizes the known 
results about the asymptotic properties and the novel results on the complexity of this law. In 
characterizing complexity, we assume diam((5) is independent of N, r and e, and we do not calculate 
how the bounds depend on r. As for the circumcenter law, we provide complete time-complexity 
results only for the case d = 1. 

Theorem 4.7 (Time and mean communication complexity of centroid law) Forr G and 
e S M+, consider the network S^d with initial conditions in Q. The following statements hold: 

(i) for d G N, the law CCccntrd achieves the e-r -deployment task '7^-r-dcpimnt/ 

(ii) for d = \ and (j) = 1, TC('7^_,.-dcplmnt, 

(Hi) for d = 1, (/) = 1 and for unidirectional communication, MCC(7^_r-dcplmnt, CCccntrd) G e(iV). . 

Proof: Fact (i) is proved in for d G {1,2} and it is clear that the same proof technique can 
be generalized to any dimension. Fact (iii) is proved in an analog OUS way to that of CCcrcmcntr ill 
Theorem 14.51 In what follows we sketch the proof of fact (ii). For d = 1, Q is a compact interval on 
M, say Q = 

We start with a brief discussion about connectivity. Note that in the r-limited Delaunay graph, 
two agents on the line that are at most at a distance r from each other are neighbors if and only if 
there are no other agents between them. Additionally we claim that, if agents i and j are neighbors 
at time instant i, then | Centroid (A' M (^)) - Centroid (^ (^))| < r. To see this, assume without loss 
of generality that xW(^) < x^^^i). Let us consider the case where the agents have neighbors on both 
sides (the other cases can be treated analogously). Let x^^{i) (respectively, x^{i)) denote the position 
of the neighbor of agent i to the left (respectively, of agent j to the right). Now, we have 

Centroid(A'[^](£)) = ^{x^^i£) + 2x^^{i) + x^\£)) , Centroid(AfW(£)) = ^(x^^ii) + 2x^^\£) + x^i\£)) . 

Therefore, | Centroid (A'W(£)) - Centroid(Af[-''] (^))| < l{\x^^{£) - x^^{£)\ + 2\x^\£) - x^^\i)\ + \x^^\£) - 

x^^{£)\) < r. This implies that agents i and j are in the same connected component of the r-limited 
Delaunay graph at time instant £ + 1. 

Next, let us consider the case that -E'r-LD(a;o) is connected. Without loss of generality, assume 
that the agents are ordered from left to right according to their identifier, that is, a;[^](0) = (a;o)i < 
••• < x^^\0) = {xo)i\f. We distinguish three cases depending on the proximity of the leftmost and 
rightmost agents 1 and N, respectively, to the boundary of the environment: (a) both agents are 
within a distance | of dQ; (b) none of the two is within a distance | of dQ; and (c) only one of 
the agents is within a distance ^ of dQ. Here is an important observation: from one time instant 
to the next one, the network configuration can fall into any of the cases described above. However, 
because of the discussion on connectivity, transitions can only occur from case (b) to either case (o) 
or (c); and from case (c) to case (a). As we show in the following, for each of these cases, the network 
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evolution under CCcontrd can be described as a discrete-time linear dynamical system which respects 
agents' ordering. 

Let us consider case (a). In this case, we have 



7(xW(^)+2xPl(^)+rE[3l(£)), ... , 



Equivalently, we can write x{ 
given by 



A. 



'+1) 



■ x[€) + ba, where the N x A^-matrix Aa and the vector ba are 





1 

4 
1 
2 



1 

4 





Note that the only equilibrium network configuration respecting the ordering of the agents is given 
by 



l))(9^ 



i £ I . 



and note that this is a |-centroidal Voronoi configuration (under the assumption of case (a)). We 
can therefore write {x{£) — x*) = Aa{x{(. — 1) — x*). Now, note that A^ = ATrid^ (i' 5) defined 
in Appendix^ Theorem IB.4r iil implies that lim^-^+oo {^i^) ~ 3^*) = Oj and that the maximum time 
required for \\x{£) — x^:\\^ < e||x(0) — x*||2 (over all initial conditions x(0) G M'^) is 0(A^^ loge^^) . It 
is not obvious, but it can be verified, that the initial condition providing the lower bound in the time 
complexity estimate does indeed have the property of respecting the agents' ordering; this fact holds 
for all three cases (a), (b) and (c). 

The case (b) can be treated in the same way. The network evolution takes now the form x{i + 1) = 
Ai, ■ x{i) + bi,, where the N x A^-matrix A^, and the vector bj, are given by 



A, 





1 
4 
1 

2 



1 
4 




4' 




L 4' 



In this case, a (non-unique) equilibrium network configuration respecting the ordering of the agents is 
of the form 



1 + N 



ir 



Note that this is a |-centroidal Voronoi configuration (under the assumption of case (b)). We can 
therefore write {x{£) — x*) = Af,{x{£ — 1) — x,,,). Now, note that Af, = ATrid^ (|, ^) as defined in 
AppendixEl As in the notation of Theorem lB.4r i) we compute Xavc = ;^l^(a^o — x^,) = jj^I'^xq. With 
this calculation, Theorem lB.4r i) implies that limf_>_|_oo {x{i) — x=k — Xavcl) = 0, and that the maximum 
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(over all initial conditions x(0) G M^) 



time required for \\x{tj — x^, — a;avcl||2 ^ ^II^(O) — x^, — Xave 
is e{N^ loge-i). 

Case (c) needs to be handled differently. Without loss of generality, assume that agent 1 is within 
distance | of dQ and agent N is not (the other case is treated analogously). Then, the network 



evolution takes now the form x(£ + 1) 
are given by 



• x{i) + be, where the x A^-matrix and the vector 



A, 



ri 

4 
1 

4 






1 

4 
1 
2 



1 
4 





Note that the only equilibrium network configuration x^^ respecting the ordering of the agents is given 
by 



1, 
+ -(2i 



l)r, 



and note that this is a |-centroidal Voronoi configuration (under the assumption of case (c)). In order 
to analyze A^, we recast the A^-dimensional discrete-time dynamical system as a 2A^-dimensional one. 
To do this, we define a 2A^-dimensional vector y by 



x^ 



',i I, and y 



[N+i] 



X 



[N-i+1] 



i £ I, 



(12) 



Now, one can see that the network evolution can be alternatively described in the variables {y^^\ ■ ■ ■ , yl^^l 
as a linear dynamical system determined by the 2N x 2N matrix ATrid^^(|, |). Using analogous argu- 
ments to the ones used before in Theorems IB . 31 and Ifj . 41 and exploiting the chain of equalities (|12jl . we 
can characterize the eigenvalues and eigenvectors of Trid2Ar_i(|, ^,\), and infer that, even for case (c), 



the maximum time required for \\x{i) — x^:\\^ < e\\x{0) — x^,||2 (over all initial conditions x(0) € 
is e{N'^ loge"^). 

In summary, for all three cases (a), (b) and (c), our calculations show that, in time O^N"^ loge~^), 
the error 2-norm satisfies the contraction inequality \\x{£) — x^\\^ < e||x(0) — rE*||2. We convert this 
inequality on 2-norms into an appropriate inequality on c«-norms as follows. Note that ||a;(0) — ||oo = 
maxjg/ |xW(0) — < (g_|_ — q-). For £ of order N'^ log ri~^ , we have: 

\\x{£) — x^Woo < \\x{i) — x^\\2 < ??||a^(0) - x^\\2 < r]VN\\x{0) — x*||oo < ?/ViV(g'+ - q~)- 



This means that e-r-deployment is achieved for ri^/N{q^ — (/_) = e, that is, in time 0{N^ logrj ^) = 
0{N^ log(iVe-i)). 

Up to here we have proved that, if the graph (/, Er-hD{xo)) is connected, then TC(7^_r-depimnti CCcentrd 
0(A''^ log(A''e~^)). If {I , Er-LD{xo)) is not connected, note that along the network evolution there can 
only be a finite number of time instants, at most — 1 where a merging of two connected components 
occurs. Therefore, the time complexity is at most 0(A^'^ log(A^e~^)). ■ 



5 Conclusions 



We have introduced a formal model for the design and analysis of coordination algorithms executed 
by networks composed of robotic agents. In our discrete-time communication, continuous-time motion 
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model, the robotic agents evolve in the physical domain in continuous-time and have the ability to 
exchange information (position and/or dynamic variables) that affect their motion at discrete-time 
instants. Under this framework, motion coordination algorithms are formalized as feedback control 
and communication laws. Drawing analogies with the classical theory on distributed algorithms, we 
have defined two measures of complexity for this formal notion of execution: the time and the mean 
communication complexity of achieving a specific task. We have defined the notion of re-scheduling of 
a control and communication law and analyzed the invariance of the proposed complexity measures 
under this operation. These concepts and results have been illustrated with various examples of 
different types of robotic networks and different classes of feedback control and communication laws. 
Finally, we have computed the proposed complexity measures for a variety of algorithms performing 
spatially-distributed tasks such as rendezvous, pursuit, and deployment. 

Throughout the paper, we have made compromises in order to arrive at a simple yet useful model 
that is able to capture a wide class of algorithms. However, it is important to acknowledge that our 
results on network, task and complexity modeling, and our results on algorithm analysis have numerous 
limitations and open up numerous avenues for future research. Let us at provide a concise list here: 

(1) modeling of asynchronous as opposed to only synchronous networks (see however [tI. I^. Eif ) : 

(2) models and analvsis of failures in the agents (arrivals/departures) and in the communication links 
(see however 0, 113, EB) 113 ) ; (3) probabilistic versions of the complexity measures, as opposed to only 
worst-case analvsis (see however jlil|): (4) quantization and delays in the communication channels 
(see however [l^ and the literature on quantized control); (5) parallel, sequential and hierarchical 
composition of control and communication laws. Additionally, our analysis results essentially consist 
of a time-complexity analysis of some basic algorithms, but many more open algorithmic questions 
remain unresolved including (6) communication complexity for omnidirectional communication models; 
(7) analysis of known algorithms for flocking and cohesion; (8) complexity analysis of tasks as opposed 
to algorithms. 
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A Basic geometric notions 

Let S C W^, d G N, be compact. The circumcenter of S, denoted by Circum(S'), is the center of the 
smallest-radius sphere in M"' enclosing S. Given an integrable function cj): S ^ M-|_, the mass of S is 
Mass(5) = Jg(p{q)dq, and the centroid of S is 



A partition of 5" is a collection of subsets of S with disjoint interiors and whose union is S. Given 
a set of N distinct points V = {Pi}iG{i,...,7v} in S, the Voronoi partition of S generated by V (with 
respect to the Euclidean norm) is the collection of sets {^i('P)}iG{i,...,Ar} defined by Vi{T') = {q G 
S \ \\q — Pi\\2 ^ Ik "Pjib ) for all S V}. We usually refer to Vi{V) as Vi. For a detailed treatment 
of Voronoi partitions we refer to 

For I = {1, . . . , N} and S C M.'^, a proximity edge map is a map of the form E: 2^^^\'^^^silxl) _ 

For r E ]R_|_, we define the r-disk proximity edge map -Er-disk '■ (K*^)^ — > '^^^^ and the r-limited Delaunay 
proximity edge map E^-hD '■ (K"^)^ 2^^^ by 



-E'r--disk(xi, . . ■,xn) = {(«, j) £ I X I\ diag(/ x /) | \\xi - Xj\\2 < r} , (13a) 
Ei^B,r{xi, ...,xn) = G / X / \ diag(/ x /) | {V^nB{xi, §)) n {VjnB{xj, §)) / 0} , (13b) 

where {Vi, . . . , Vn} is the Voronoi partition of M*^ generated by {xi, . . . ,xj\i}. Illustrations of these 
concepts are given in Fig. d For d= 1, the r-limited Delaunay proximity edge map has the following 
intuitive characterization: two points are neighbors if and only if they are within distance r and no 
other point is between them. 

As proved in 0], the r-limited Delaunay graph and the r-disk graph have the same connected 
components. Additionally, the r-limited Delaunay graph is "computable" on the r-disk graph in 
the following sense: any node in the network can compute the set of its neighbors in the r-limited 
Delaunay graph if it is given the set of its neighbors in the r-disk graph. This implies that any control 
and communication law for a network with communication graph Er-hD can be implemented on a 
analogous network with communication graph -Br-disk- 
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Figure 7: The r-disk and r-limited Delaunay graphs in 

B Tridiagonal Toeplitz and circulant matrices 

We briefly review some basic facts about certain classes of Toeplitz matrices, see For N >2 and 
a,b,c € M, define the N x N Toeplitz matrices TridAr(a, b, c) and CircAr(a, b, c) by 



TridAr(a, 5, c) 



b c 
a b c 



a b c 
a b 



, and CircAr(a, 5, c) = TridAr(a, 5, c) + 



... 

... 



c 



a 









The matrices TridAr and CircAr are tridiagonal and circulant, respectively. The two matrices only differ 
in their (1, N) and {N, 1) entries. Note our convention that 6*2(0, b,c) = ^ ~J)'^ 

Lemma B.l (Eigenvalues of tridiagonal Toeplitz and circulant matrices) For N > 2 and 

a,b,c £ M, the following statements hold: 



(i) for ac 7^ 0, the eigenvalues and eigenvectors 0/ TridAr (a, 6, c) are 



6 + 2c^/- cos ! 

c V ^ + 1 



, {!,... ,iV}; 



(ii) the eigenvalues and eigenvectors of Circj\f {a, b,c) are (note oo = exp(^^^^^^)j 



/i27r\ , / f27r\ 

6+ (a + c)cos ( — I + V-l(c- a)sm I — I 



1 

{N-l)i 



{!,..., iV}. 



Proof: Both facts are discussed, for example, in |3ll . Example 7.2.5 and Exercise 7.2.20]. Fact (ii) 
requires some straightforward algebraic manipulations. ■ 

Remarks B.2 (i) The set of eigenvalues of TridAr(a, ^, c) is contained in the real interval [6 — 
2-^/ac, 6 + 2-^/ac], if ac > 0, and in the interval in the complex plane [b — 2\/—ly^\ac\, b + 
2y/^J\^\], if ac < 0. 
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(ii) The set of eigenvalues of CircAr(a, 6, c) is contained in the ehipse on the complex plane with 
center b and horizontal axis 2\a + c\ and vertical axis 2|c — a\. 



(iii) Recall from [3l| that (1) a square matrix is normal if it has a complete orthonormal set of 
eigenvectors, (2) circulant matrices and real-symmetric matrices are normal, and (3) if a normal 
matrix has eigenvalues {Ai, . . . , A„}, then its singular values are {|Ai|, . . . , |Ari|}. • 

We can now state the main result of this section. 

Theorem B.3 (Tridiagonal Toeplitz and circulant dynamical systems) Let N > 2, e g]0, 1[, 
and a, 6, c G M. Let x: 'Nq ^ and y : No ^ be solutions to 

x{£ + 1) = TTidiy{a,b,c) x{£), and y{i + I) = Cii[ciy{a,b,c) y{£), 

with initial conditions x{0) = xq and y(0) = yo, respectively. The following statements hold: 



(i) if a = c ^ and \b\ + 2\a\ = 1, then lim^ 



+ + 00 



x{i) = 0, and the maximum time required for 



^ ^ll^^olb (over all initial conditions xq G 



IS 



(ii) if a ^ 0, c = and < |6| < 1, then lim^^+oo = 0, and the maximum time required for 
\\x{i)\\2 < ell^^olb (over all initial conditions xq E M^J is 0(A^log + loge"-^); 



N 



(Hi) if a > 0, c > 0, b > 0, and a + 6 + c = 1, then limg^^^y{£) = yavcl, where yav( 

and the maximum time required for \\y{£) — yavcllb ^ ^||yo ~ ?/avcl||2 (over all initial conditions 
yoeR^) is e{N^ log e~^). • 

Proof: Let us prove fact (i). We start by bounding from above the eigenvalue with largest 
absolute value, that is, the largest singular value, of TridAr(a, 6, a): 



max 

ie{i,...,N} 



b + 2a cos 



ITT 



N +1 



< 161 + 2\a\ max 

i&{l,...,N} 



cos 



N +1 



< \b\ + 2\a\ cos 



TT 



iV+ 1 



Because cos(-^yq;Y) < 1 for any N > 2, the matrix TridAr(a, 6, a) is stable. Additionally, for ^ > 0, we 
bound from above the magnitude of the curve x as 



TridAr(a, b, a) xo\\2 < ( |6| + 2\a\ cos 



TT 



iV + 1 



|2;o||2- 



In order to have ||a;(£)||2 < ella^olb) it is sufficient that £log + 2\a\ cos (^^jy-qTyj j < loge, that is 

loge"^ 



> 



log(|6| + 2|a|cos (jv^ 



(14) 



To show the upper bound, note that as t — > we have 

1 



1 



+ 0(1). 



log(l - 2|a|(l - cost)) \a\t^ 
Now, assume without loss of generality that ab > and consider the eigenvalue b + 2a cos{ ^^_^-^ ) of 



TridAr(o, 6, a). Note that |6+2a cos( 



b + 2ocos( 



Ntv 
N+l 



N+1' 



|6|+2|a| cos( 



(If ab < 0, then consider the eigenvalue 



I.) For N > 2, define the unit-length vector 



VAT 



N + l 



sm 



sm 



N+l 

Nn 
N+l. 



oN 



(15) 
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and note that, by Lemma lB.ir i). vjv is an eigenvector of TridAr(a, b, a) with eigenvalue b + 2a cos{ j^^^ ). 

The trajectory x with initial condition v^r satisfies ||x(£)||2 = + 2|a| cos ^■jyq:^-^^ l|vAr||2 and, 
therefore, it will enter 5(0, e||vAr||2) only when i satisfies (|14|) . This completes the proof of fact (i). 
Next we consider statement (ii). Clearly, TridAr(a, 6, 0) is stable. For £ > 0, we compute 



N~l 



Trid^(a, b, OY = b' (/^ + ^ TVid^(l, 0, 0)) = b' (^)' TVid;v(l, 0, 0)^' 
because of the nilpotency of Tridjv(l, 0, 0). Now we can bound from above the magnitude of the curve 



X as 



7V-1 

2 = ||Trid;v(a,?^,0)'3;o||2 < l^l' E ^ 



j)\\b 



a\'J 



Tridjv(l, 0,0)^X0 1 



<e'^/^£^-i|6n|xo||2. 



Here we used || TridAr(l, 0, Oyxolb < ||xo||2 and maxj^j^jyi I J e {0, • • • , iV - 1}} < Therefore, 
in order to have ||x(^)||2 < e||a;o||2, it suffices that log(e°/'') + {N — 1) log£ + ^log |6| < loge, that is 



N-l f -loge 

log£ > 



log|6| 



log|6| 



A sufficient condition for ^ — a log i > (3, for a, /3 > 0, is that £ > 2/3 + 2a max{l, log a}. For, if £ > 2a, 
then log£ is bounded from above by the line i/2a + log a. Furthermore, the line i/2a + log a is a 
lower bound for the line {i — f3)/a if £ > 2/3 + 2a log a. In summary, it is true that ||a;(£)||2 < e||a:;(0)||2 
whenever 



> 2-^ 



log|5| 



+ 2- 



A^- 1 
log|6| 



max < 1, log ■ 



N- 1 
-log|6| 



This completes the proof of the upper bound, that is, fact (ii). 

The proof of fact (iii) is similar to that of fact (i). We analyze the singular values of CircAr(a, b, c). 
It is clear that the eigenvalue corresponding to i = is equal to 1; this is the largest singular value 
of CircAr(a, 6, c) and the corresponding eigenvector is 1. In the orthogonal decomposition induced by 
the eigenvectors of CircAr(a, b, c), the vector i/q has a component ?/ave along the eigenvector 1. We now 
compute the second largest singular value: 



max 

i&{l,...,N~l} 



i27r\ 



6 + (a + c) cos — + V^(c - a) sin — 



A / 



i2TT\ 



A 7 



,27r 



1 - (a + c) 1 - cos ( — ) + V^(c - a) sin — 



A 



27r\ 



A J 



Here | • | is the norm in C. Because of the assumptions on a, b, c, the second largest singular value is 
strictly less than 1. For £ > 0, we bound the distance of the curve y{£) from yavel as 



\y{^) - yavcllb = II CircAr(a, b, cfy^ - yavcllb 

27r^ 



CircAr(a, c)^(yo - l/avcl) lb 



< 



1 - (« + c) (l - COS (-^)) + \/^(c - a) sin i 



2tt\ 



\yo - yavellb- 
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This proves that hm£_^+oo y(^) = ^avcl- Also, for a = a + c, P = c — a and as t — > 0, we have 

1 2 



log ( (1 - q(1 - cost)) +(3^sm'^t^ ^ ' 



+ 0(1). 



Here 0^ < a because a,c g]0, 1[. 



Now, consider the eigenvalues \n = h + {a + c) cos y-^) + ^/^^{c — a) sin (-^) and A^r = h + 

of CircAr(a, 6, c), and its associated eigenvectors (cf. 



(a + c) cos I -^^ — — I 
Lemma IB.ll iil ) 



+ V— l(c — a) sin 



(jV-l)27r 
N 



VAT 



W 



1 

If 



w 



1 

TV-l 



G C 



TV 



(16) 



Note that the vector v^v + v^r belongs to M^. Moreover, its component yave along the eigenvector 
1 is 0. The trajectory y with initial condition va? + v^r satisfies ||y(^)||2 = HA^rVTv + AArVAr||2 = 
|AAr|^||vAr + V7v||2 and, therefore, it will enter B{Q,e\\wN + v^vlb) only when 

loge"^ 



log 



1 - (a + c) ( 1 - cos (^) ) + ^f=l{c - a) sin (f 



This completes the proof of fact (iii). ■ 
Next, we extend these results to another interesting set of matrices. For N > 2 and a, 6 £ M, define 
the N X N matrices ATrid^(a, 6) and ATrid]^(a,fe) by 



If we define 













a 


... 




0" 





















... 

















ATrid 




Tridjv(a, 6, a) ± 





























































a 








1 


1 





... 


" 




1 


1 













1 


-1 


1 


... 







-1 


1 


1 










1 





-1 


1 ... 





, and P- 




1 





1 


1 







1 







-1 


1 




(_l)TV-2 










1 


1 


1 










-1 

















1 



then the following similarity transforms are satisfied: 

'b + 2a 

TridAr_i(a, 6, a) 

'b-2a 

TridAr_i(a, 6, a) 



ATrid^(a,6) = P^ 
ATrid]^(a,6) = P_ 



p-i 
P~\ 



(17) 



To analyze the convergence properties of the dynamical systems determined by ATrid^(a, 6) and 

ATrid^(a, 6), we recall that 1 = (1, . . . , 1) G M"^, and we define 1_ = (1, -1, 1, . . . , (-1)^-^^ (-1)^-^) e 
Tn>N 
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Theorem B.4 Let N >2, e g]0, 1[, and a,b€R with a ^ and \b\ +2\a\ = 1. Letx: No and 
2; : No — > be solutions to 

x{i + 1) = ATrid^(a, b) x{i), and z{£ + 1) = ATrid^(a, b) z{l), 

with initial conditions x{0) = xq and z(0) = zq, respectively. The following statements hold: 

(i) lim£_>_|_oo {x{t) — j;ave(^)l) = 0, where Xavel^) = (;^1^2;o)(6 + 2a)^, and the maximum time 
required for \\x{£) — Xave(^)l||2 < ^lla^o — 3;ave(0)l||2 (over all initial conditions xq G M^J is 
9(^2 loge-i); 

(a) {z{£) — Zavc{^)'^~) = 0; where z^vei^) = {jj''^'^zo){b — 2aY , and the maximum time 

required for \\z{i) — 5;ave(^)l- II2 < ^W^o — -2ave(0)l- II2 (over all initial conditions zq G M^J is 
e(iV2 loge-i). 

Proof: We prove fact (i) and remark that the proof of fact (ii) is analogous. Consider the change 
of coordinates 



x{e) = P4 



x'^^^{£)l + P^ 







where Xave(^) G ^ and y{£) G ^. A quick calculation shows that ^aveC^) 
similarity transformation described in equation H17() implies 



jfl^x{£), and the 



y(£ + 1) = Trid,v-i(a, b, a) y{£), and x',,,i£ + 1) = (b + 2a)x'^,,i£). 



Therefore, x^vc 



x(£+l) 



x'avc- It is also clear that 





y{£+l) 





TridAr_i(a, a) 



P+l)(x(£)-XavcWl). 



Consider the matrix in parenthesis determining the trajectory £ 1— > {x{£) — a;ave(-^)l)- This matrix is 
symmetric, its singular values are and the singular values of Trid7v_i(a, 6, a), and its eigenvectors 
are 1 and the eigenvectors of TridAr_i(a, 6, a) (padded with an extra zero). These facts are sufficient 
to duplicate, step by step, the proof of fact (i) in Theorem IB. 31 Therefore, the trajectory £ 1— > 
{x{£) — Xave(-^)l) satisfies the same properties as those stated in Theorem IB.Sr i). ■ 
We conclude this section with some useful bounds. 

Lemma B.5 Assume x G R^, y G R^~-^ and z G M^"-*^ jointly satisfy 



Then i||x||2 < ||y||2 < (iV - l)||x||2 and ^\\x\\2 < \\z\\2 < {N - l)\\x\\2. 
Proof: The proof is based on the coordinate expressions: 

xi = yi, X2 = y2- yi, ■■■ xn-i = yN-i - yN-2, xn = -vn-i, 

yi=Xi, y2 = X2+Xi, 2/3 = X3 + X2 + Xi, ... yN-1 = XN-1 H h Xi, 

and 

Xi = Zi, X2 = Z2 + Zi, ... XAr_i = Zj^-l + ZN-2, Xn = Z^-l, 
Zi= Xi, Z2 = X2 - Xi, Z3 = X3 - X2 + Xi, ... ZN-1 = Xn~1 H h (-1)^"^ 



Xi. 
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